Summary. It is explained how field-theoretic methods and the dynamic renor- 
malisation group (RG) can be applied to study the universal scaling properties of 
systems that either undergo a continuous phase transition or display generic scale in- 
variance, both near and far from thermal equilibrium. Part 1 introduces the response 
functional field theory representation of (nonlinear) Langevin equations. The RG is 
employed to compute the scaling exponents for several universality classes governing 
the critical dynamics near second-order phase transitions in equilibrium. The effects 
of reversible mode-coupling terms, quenching from random initial conditions to the 
critical point, and violating the detailed balance constraints are briefly discussed. It 
is shown how the same formalism can be applied to nonequilibrium systems such as 
driven diffusive lattice gases. Part 2 describes how the master equation for stochas- 
tic particle reaction processes can be mapped onto a field theory action. The RG is 
then used to analyse simple diffusion-limited annihilation reactions as well as generic 
continuous transitions from active to inactive, absorbing states, which are charac- 
terised by the power laws of (critical) directed percolation. Certain other important 
universality classes are mentioned, and some open issues are listed. 
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1 Critical Dynamics 

Field-theoretic tools and the renormalisation group (RG) method have had 
a tremendous impact in our understanding of the universal power laws that 
emerge near equilibrium critical points (see, e.g., Refs. El E El El ) j in- 

cluding the associated dynamic critical phenomena jJlEl. Our goal here is to 
similarly describe the scaling properties of systems driven far from thermal 
equilibrium, which either undergo a continuous nonequilibrium phase transi- 
tion or display generic scale invariance. We are then confronted with cap- 
turing the (stochastic) dynamics of the long-wavelength modes of the 'slow' 
degrees of freedom, namely the order parameter for the transition, any con- 
served quantities, and perhaps additional relevant variables. In these lecture 
notes, I aim to briefly describe how a representation in terms of a field theory 
action can be obtained for (1) general nonlinear Langevin stochastic differen- 
tial equations [HI El; and (2) for master equations governing classical particle 
reaction-diffusion systems ^1 ^2 ^1 ■ I will then demonstrate how the dy- 
namic (perturbative) RG can be employed to derive the asymptotic scaling 
laws in stochastic dynamical systems; to infer the upper critical dimension 
dc (for dimensions d < d^, fluctuations strongly affect the universal scaling 
properties); and to systematically compute the critical exponents as well as to 
determine further universal properties in various intriguing dynamical model 
systems both near and far from equilibrium. (For considerably more details, 
especially on the more technical aspects, the reader is referred to Ref. 13j.) 

1.1 Continuous phase transitions and critical slowing down 

The vicinity of a critical point is characterised by strong correlations and large 
fluctuations. The system under investigation is tlien behaving in a highly coop- 
erative manner, and as a consequence, the standard approximative methods of 
statistical mechanics, namely perturbation or cluster expansions that assume 
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either weak interactions or short-range correlations, fail. Upon approaching 
an equilibrium continuous (second-order) phase transition, i.e., for |r| ^ 1, 
where r = (T — T^) /T^ measures the deviation from the critical temperature 
Tc, the thermal fluctuations of the order parameter S{x) (which characterises 
the different thermodynamic phases, usually chosen such that the thermal 
average (S) = vanishes in the high-temperature 'disordered' phase) are, in 
the thermodynamic limit, governed by a diverging length scale 

e(r) ^ \rr ■ (1) 

Here, we have defined the correlation length via the typically exponential 
decay of the static cumulant or connected two-point correlation function 
C{x) = {S{x)S{0)) — (S)^ ~ e^l'^l/?, and v denotes the correlation length 
critical exponent. As T —+ Tc, ^ —^oo, which entails the absence of any charac- 
teristic length scale for the order parameter fluctuations at criticality. Hence 
we expect the critical correlations to follow a power law C{x) ^ 
in d dimensions, which defines the Fisher exponent rj. The following scaling 
ansatz generalises this power law to T ^ Tc, but still in the vicinity of the 
critical point, ^ 

CiT,x) = \x\-^^~^+^^C±ix/0 , (2) 

with two distinct regular scaling functions (7+ (y) for T > and C_ (y) for 
T < Tc, respectively. For its Fourier transform C(r, q) — J d'^x e~"^'^ C(t, a;), 
one obtains the corresponding scaling form 

C{T,q) = \q\-^+'^C±{qO , (3) 

with new scaling functions C±{p) = jpp-'' Jd''ye~'P y |y|-('i-2+'?) C±{y). 

As we will see in Section 11.51 there are only two independent static crit- 
ical exponents. Consequently, it must be possible to use the static scaling 
hypothesis (jUl or along with the definition QJ, to express the expo- 
nents describing the thermodynamic singularities near a second-order phase 
transition in terms of v and rj through scaling laws. For example, the or- 
der parameter in the low-temperature phase (t < 0) is expected to grow 
as (S) ^ {—t}'^. Let us consider Eq. ^ in the limit \x\ oo. In or- 
der for the \x\ dependence to cancel, C±{y) oc ly]"^^'^^^ for large \y\, and 
therefore C{t,\x\ -> cx)) - ^-{d-2+n) _ |^|i'(<i-2+,,)^ Qn the other hand, 

C(t, IccI -> oo) -(S')2 (-r)2^ for T < T^; thus we identify the order 

parameter critical exponent through the hyper scaling relation 

P^'^id-2 + r^) . (4) 

Let us next consider the isothermal static susceptibility Xn which according to 
the equilibrium fluctuation-response theorem is given in terms of the spatial 
integral of the correlation function C(t, x): Xt{t) — (^b^)"^ limq_^o C{t, q). 
But C±{p) ^ as p ^ to ensure nonsingular behaviour, whence 
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XriT) ^ (.^ ^ \t\ ^'^'^ ^\ and upon defining the associated thermodynamic 
critical exponent 7 via XtIt) ^ I'^'l^^i we obtain the scaling relation 

7 = ^(2-r;). (5) 

The scaling laws Q , ® as well as scaling relations such as Q and ||SJ| can 
be put on solid foundations by means of the RG procedure, based on an ef- 
fective long- wavelength Hamiltonian Ti.[S], a functional of S{x), that captures 
the essential physics of the problem, namely the relevant symmetries in order 
parameter and real space, and the existence of a continuous phase transition. 
The probability of finding a configuration S{x) at given temperature T is then 
given by the canonical distribution 

Pcq[5] ocexp(-7^[5]/fcBr) . (6) 

For example, the mathematical description of the critical phenomena for an 
0(n)-symmetric order parameter field S'^^x), with vector index a = 1, . . . , rt, 
is based on the Landau-Ginzhurg-Wilson functional ^ 12 El El ISl El 



n[s] = j A'^x 



'-[S-{xr + \[VS-{x)f 



--[S%x)YY.^SP{x)f-h-{x)S-{x) 



(7) 



where h°'{x) is the external field thermodynamically conjugate to S°'{x), 
u > denotes the strength of the nonlinearity that drives the phase transfor- 
mation, and r is the control parameter for the transition, i.e., r cx T — T^, 
where is the (mean-field) critical temperature. Spatial variations of the 
order parameter are energetically suppressed by the term ~ [VS'"(a;)]^, and 
the corresponding positive coefRcient has been absorbed into the fields S". 

We shall, however, not pursue the static theory further here, but instead 
proceed to a full dynamical description in terms of nonlinear Langevin equa- 
tions [3 El ■ We will then formulate the RG within this dynamic framework, 
and therein demonstrate the emergence of scaling laws and the computation 
of critical exponents in a systematic perturbative expansion with respect to 
the deviation e = d — dc from the upper critical dimension. 

In order to construct the desired effective stochastic dynamics near a crit- 
ical point, we recall that correlated region of size ^ become quite large in the 
vicinity of the transition. Since the associated relaxation times for such clus- 
ters should grow with their extent, one would expect the characteristic time 
scale for the relaxation of the order parameter fluctuations to increase as well 
as T —^ Tc, namely 

ic(r)^^(rr^|r|— , (8) 

which introduces the dynamic critical exponent z that encodes the critical 
slowing down at the phase transition; usually z > 1. Since the typical relax- 
ation rates therefore scale as uJc{t) = l/tdr) ~ \t\^^ , we may utilise the static 
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scaling variable p — to generalise the crucial observation ||SJ) and formulate 
a dynamic scaling hypothesis for the wavevector-dependent dispersion relation 
of the order parameter fluctuations |14l I15| . 

u;c{T,q)^\q\'LO±{qO ■ (9) 

We can then proceed to write down dynamical scaling laws by simply 
postulating the additional scaling variables s = t/tdr) or uj/ujc{T,q). For 
example, as an immediate consequence we find for the time-dependent mean 
order parameter 

{S{T,t))^\T\l'S{t/Q , (10) 

with S{s ~f oo) — const., but S{s) ~ g-Z^/^i^ as s — + in order for the r 
dependence to disappear. At the critical point {t — 0), this yields the power- 
law decay {S{t)) ^ i^", with 

. = A._L(,_2 + ,). (11) 

Similarly, the scaling law for the dynamic order parameter susceptibility ( re- 
sponse function ) becomes 

x(T,q,c.) = |q|-2+''x±(9C,c.r) , (12) 

which constitutes the dynamical generalisation of Eq. J^Jl, for x(t, q,0) = 
{k'BT)~^ C {t , q) . Upon applying the fluctuation-dissipation theorem, valid in 
thermal equilibrium, we therefrom obtain the dynamic correlation function 

C{T,q,u;) = ^ lmx{r,q,Lo) = \qr''+^ C±{q^,Lun , (13) 

UJ 

and for its Fourier transform in real space and time, 

C{r,x,t)^ I -0-^ l^e'^''----'^ C{r,q,u;) = \x\-^''-'+^'^ C±ix/U/n , 

(14) 

which reduces to the static limit Q if we set t — 0. 

The critical slowing down of the order parameter fluctuations near the crit- 
ical point provides us with a natural separation of time scales. Assuming (for 
now) that there are no other conserved variables in the system, which would 
constitute additional slow modes, we may thus resort to a coarse-grained long- 
wavelength and long-time description, focusing merely on the order parameter 
kinetics, while subsuming all other 'fast' degrees of freedom in random 'noise' 
terms. This leads us to a mesoscopic Langevin equation for the slow variables 
S°'(x,t) of the form 



(15) 
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In the simplest case, the systematic force terms here just represent purely 
relaxational dynamics towards the equilibrium configuration , 

where D represents the relaxation coefficient, and 7i[5] is again the effective 
Hamiltonian that governs the phase transition, e.g. given by Eq. iQ. For the 
stochastic forces we may assume the most convenient form, and take them to 
simply represent Gaussian white noise with zero mean, (^"(a;,i)) = 0, but 
with their second moment in thermal equilibrium fixed by Einstein's relation 

{C{x,t)C'^{x',t')) = 2kBTDS{x-x')S{t-t')S°''^ . (17) 

As can be verified by means of the associated Fokkcr-Planck equation for 
the time-dependent probability distribution V [S, t] , Eq. (|17|1 guarantees that 
eventually 7'[S', i ^ oo] Vcq[S], the canonical distribution The stochas- 
tic differential equation (|15|l . with H16|l . the Hamiltonian and the noise 
correlator (|17|l , define the relaxational model A (according to the classification 
in Ref. [7]) for a nonconserved 0(n)-symmetric order parameter. 

If, however, the order parameter is conserved, we have to consider the as- 
sociated continuity equation dt S°' + 'V ■ J" = 0, where typically the conserved 
current is given by a gradient of the field S": J" = — D 'V 3°" + . . .; as a con- 
sequence, the order parameter fluctuations will relax diffusively with diffusion 
coefficient D. The ensuing model B 16, 7^ for the relaxational critical dynam- 
ics of a conserved order parameter can be obtained by replacing D —D 
in Eqs. 116(1 and (|17l) . In fact, we will henceforth treat both models A and B 
simultaneously by setting D D (iV)'', where a — and a — 2 respectively 
represent the nonconserved and conserved cases. Explicitly, we thus obtain 

ds'^jx^t) snjs] 

+D{iVyh''{x,t) + C{x,t), (18) 



with 



{C{x,t)(:'^{x',t')) =2kBTD{iVY5[x-x')5[t-t')5°'^ . (19) 



Notice already that the presence or absence of a conservation law for the order 
parameter implies different dynamics for systems described by identical static 
behaviour. Before proceeding with the analysis of the relaxational models, we 
remark that in general there may exist additional reversible contributions to 
the systematic forces F°'[S], see Sec. 11.61 and / or dynamical mode-couphngs 
to additional conserved, slow fields, which effect further splitting into several 
distinct dynamic universality classes [f)l 171 IT^. 
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Let us now evaluate the dynamic response and correlation functions in the 
Gaussian (mean-field) approximation in the high-temperature phase. To this 
end, we set u = and thus discard the nonlinear terms in the Hamiltonian 
(Q as well as in Eq. H18|) . The ensuing Langevin equation becomes linear in 
the fields S", and is therefore readily solved by means of Fourier transforms. 
Straightforward algebra and regrouping some terms yields 

[-iu; + Dq- {r + q')] lo) = Dq'^ h"{q, u:) + ^ (g, u:) . (20) 

With (C"(q,a;)) = 0, this gives immediately 

= i?q'^Go(g,c.)5"^ (21) 

with the response propagator 

Go{q, Lo) = [-ice; + Dq'^ (r + q")] . (22) 

As is readily established by means of the residue theorem, its Fourier back- 
transform in time obeys causality^ 

Go(q,i) -©(t)e-^''" . (23) 

Setting h" = 0, and with the noise correlator (|19|l in Fourier space 

(r (q, Lo) C^{q', c.')) = 2fcBT Dq^ {2Trf+^6{q + q') 5{u + lo') 5^^ , (24) 

we obtain the Gaussian dynamic correlation fimction (^S°'{q, uj) S^{q' , w'))o ~ 
Co(q, lo) (27r)''+iJ(q + q') 5{lo + lo'), where 

= c.^ + [^g^(?+g^)P - 2^^^^'^" \Go{q,^)\' ■ (25) 

The fluctuation-dissipation theorem H13|l is of course satisfied; moreover, as 
function of wavevector and time, 

Co(q,i) = -^e-^''"('-+''')l*l . (26) 
r + q^ 

In the Gaussian approximation, away from criticality (r > 0, q 7^ 0) the tem- 
poral correlations for models A and B decay exponentially, with the relaxation 
rate LOc{r,q) = I?q^"'"°(l -I- r/q^). Upon comparison with the dynamic scal- 
ing hypothesis Q, we infer the mean-field scaling exponents vq = 1/2 and 
2:0 = 2 + a. At the critical point, a nonconserved order parameter relaxes diffu- 
sively (zq = 2) in this approximation, whereas the conserved order parameter 
kinetics becomes even slower, namely subdiffusive with zq = 4. Finally, invok- 
ing Eqs. (O, (inil, CH), or simply the static limit Co(q, 0) = k^T/lr + q^), 
we find r/o = for the Gaussian model. 



Q/3 

Xo 



djS'^jq^Lo)) 
dhP{q,uj) 
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The full nonlinear Langevin equation H18|) cannot be solved exactly. Yet 
a perturbation expansion with respect to the coupling u may be set up in a 
slightly cumbersome, but straightforward manner by direct iteration of the 
equations of motion I17j . More elegantly, one may utilise a path-integral 
representation of the Langevin stochastic process , which allows the 

application of all the standard tools from statistical and quantum field theory 
0121 El EI El m 7 and has the additional advantage of rendering symmetries in 
the problem more explicit (HUHIEI- 



1.2 Field theory representation of Langevin equations 



Our starting point is a set of coupled Langevin equations of the form (|15|l for 
mesoscopic, coarse-grained stochastic variables S'^{x,t). For the stochastic 
forces, we make the simplest possible assumption of Gaussian white noise, 

{C{x,t))=0, {C{x,t)C'^{x',t')) ^2L°'6{x~x')6{t~t')5°'^ , (27) 

where L" may represent a differential operator (such as the Laplacian for 
conserved fields) , and even a functional of 3°" . In the time interval < t < tf, 
the moments H27|l are encoded in the probability distribution 



>V[C] cx exp 



i / d'x f'dt ^ C (x, t) [{L'-r'C (x, t)] 



(28) 



If we now switch variables from the stochastic noise to the fields S"' by 
means of the equations of motion (|15|) , we obtain 

>V[C] I?[C] = 'PiS] T)[S] oc e-^[^l V[S] , (29) 

with the statistical weight determined by the Onsager-Machlup functional |2] 



Q[S]^\jA'xjdtY^{^-^^F-[S] 



(30) 



Note that the Jacobian for the nonlinear variable transformation {C"} ~^ 
{S*"} has been omitted here. In fact, the above procedure is properly defined 
through appropriately discretising time. If a forward (Ito) discretisation is 
applied, then indeed the associated functional determinant is a mere constant 
that can be absorbed in the functional measure. The functional H30|) already 
represents a desired field theory action. Since the probability distribution for 
the stochastic forces should be normahsed, /l?[C] W[Q = 1, the associated 
'partition function' is unity, and carries no physical information (as opposed 
to static statistical field theory, where it determines the free energy and hence 
the entire thermodynamics). The Onsager-Machlup representation is however 
plagued by technical problems: Eq. contains (i")~^, which for conserved 
variables entails the inverse Laplacian operator, i.e., a Green function in real 
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space or the singular factor in Fourier space; moreover the nonhnearities 
in F°' [S] appear quadratically. Hence it is desirable to linearise the action H3U|) 
by means of a Hubbard-Stratonovich transformation P|. 

We shall follow an alternative, more general route that completely avoids 
the appearance of the inverse operators (L")~^ in intermediate steps. Our 
goal is to average over noise 'histories' for observables A[S] that need to be 
expressible in terms of the stochastic fields S": {A[S])i cx JV[C] A[S{C)] W[C]. 
For this purpose, we employ the identity 



a (x,t) 

V[iS] jV[S] exp 



dS°'(x,t) 
di 



F^[S]{x,t)-CM 
dS°' 



I V[iS] I V[S] exp -Jd'xJ dt ^ 5" - - C" 



, (31) 



where the first line constitutes a rather involved representation of the unity (in 
a somewhat symbolic notation; again proper discretisation should be invoked 
here), and the second line utilises the Fourier representation of the (functional) 
delta distribution by means of the purely imaginary auxiliary fields S (and 
factors 2tt have been absorbed in its functional measure). 

Inserting and the probability distribution into the desired stochas- 
tic noise average, we arrive at 



{A[S])(; CX IV[iS] / Pis'] exp 
'P[C]exp( - f d'^x 



A[S] 
(32) 



We may now evaluate the Gaussian integrals over the noise C", which yields 



{A[S])c = V[S]A[S]V[S] , V[S] 



V[iS] < 



-A[S,S] 



(33) 



with the statistical weight now governed by the Janssen-De Dominicis 're- 
sponse 'functional [THIIBI I^ 



A[S,S] = j d'^xj^ ^dt^ 



ga j^a ga 



(34) 



Once again, we have omitted the functional determinant from the variable 

change ^ {5*"}, and normaHsation impHes /^[iS] /X»[S'] e--^[^"51 = 1. 
The first term in the action (|34|l encodes the temporal evolution according 
to the systematic terms in the Langevin equations (|15|l . whereas the second 
term specifies the noise correlations (|27ll . Since the auxiliary variables 5*", 
often termed Martin-Siggia-Rose response fields [201, appear only quadrati- 
cally here, they may be eliminated via completing the squares and Gaussian 
integrations; thereby one recovers the Onsager-Machlup functional (|30|l . 



Field Theory Approaches to Nonequilibrium Dynamics 



11 



The Janssen-De Dominicis functional (|34|) takes the form of a (d + 1)- 
dimensional statistical field theory with two independent sets of fields S" and 
S" . We may thus bring the established machinery of statistical and quantum 
field theory m^lElEElini to bear here; it should however be noted that the 
response functional formalism for stochastic Langevin dynamics incorporates 
causality in a nontrivial manner, which leads to important distinctions [H'. 

Let us specify the Janssen-De Dominicis functional for the purely relax- 
ational models A and B see Eqs. (fH^ and (f11?| . splitting it into the 

Gaussian and anharmonic parts A = .4o + -4int |9j , which read 



Ao[S,S] = j cV^x j dt^yS"" 



Aint [S,S]=D^ J d'^x JdtJ^S" (i V)"^ 5^ Sf^ . 



(35) 
(36) 



a, 13 



Since we are interested in the vicinity of the critical point T ^ Tc, we have 
absorbed the constant kBTc into the fields. The prescription (|33|l tells us how 
to compute time-dependent correlation functions (^S°'{x,t) S^{x' ,t')y Using 
Eq. (|35|l . the dynamic order parameter susceptibility follows from 



S{S''{x,t)) 



5hf^{x',t') 



^D(s^{x,t){iVrsP{x',t')) ; 



h=Q 



(37) 

for the simple relaxational models (only), the response function is just given by 
a correlator that involves an auxiliary variable, which explains why the S°' are 
referred to as 'response' fields. In equilibrium, one may employ the Onsager- 
Machlup functional H3()(l to derive the fluctuation-dissipation theorem [HI 



d 



X'^^ix -x\t- t') = 0{t - t') ^ {S'^ix, t) S^{x', t')) 



(38) 



which is equivalent to Eq. in Fourier space. 

In order to access arbitrary correlators, we define the generating functional 



z[j,j] = ^exp J d^x Idtj2 {r 5" +r 



(39) 



wherefrom the correlation functions follow via functional derivatives. 



(40) 



j=o=j 



and the cumulants or connected correlation functions via 



12 Uwe Claus Tauber 



In the harmonic approximation, setting u — 0, Z[j,j] can be evaluated ex- 
plicitly (most directly in Fourier space) by means of Gaussian integration 
IHIEI; one thereby recovers (with k^T — 1) the Gaussian response propaga- 
tor (|22|l and two-point correlation function (|25|l . Moreover, as a consequence 

of causality, (s"iq,uj) S^^(q', w'))^ = 0. 



(41) 



1.3 Outline of dynamic perturbation theory 

Since we cannot evaluate correlation functions with the nonlinear action H36|l 
exactly, we resort to a perturbational treatment, assuming, for the time being, 
a small coupling strength u. The perturbation expansion with respect to u is 
constructed by rewriting the desired correlation functions in terms of averages 
with respect to the Gaussian action (|35|l . henceforth indicated with index '0', 
and then expanding the exponential of —Amt, 



Yi s*"' 



1 

°° 1 i \ 

\ . (42) 



The remaining Gaussian averages, a series of polynomials in the fields S"" 
and 5'", can be evaluated by means of Wick's theorem, here an immediate 
consequence of the Gaussian statistical weight, which states that all such av- 
erages can be written as a sum over all possible factorisations into Gaussian 
two-point functions (^S" S^'^^, i.e., essentially the response propagator Gq, 
Eq. ((221), and {S" S^^)^, the Gaussian correlation function Co, Eq. ((23 • Re- 
call that the denominator in Eq. H42|l is exactly unity as a consequence of 
normalisation; alternatively, this result follows from causality in conjunction 
with our forward descretisation prescription, which implies that we should 
identify 0(0) = 0. (We remark that had we chosen another temporal dis- 
cretisation rule, any apparent contributions from the denominator would be 
precisely cancelled by the in this case nonvanishing functional Jacobian from 
the variable transformation {C"} {S"}-) At any rate, our stochastic field 
theory contains no 'vacuum ' contributions. 

The many terms in the perturbation expansion H42|) are most lucidly or- 
ganised in a graphical representation, using Feynman diagrams with the basic 
elements depicted in Fig. ^ We represent the response propagator by a 
directed line (here conventionally from right to left) , which encodes its causal 
nature; the noise by a two-point 'source' vertex, and the anharmonic term 
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(a) 



a 



P - ia+ D q''(r + q'^) 



(b) 



(c) 



, a 



P = - D?"-^ 
o 



Fig. 1. Elements of dynamic perturbation theory for the 0(n)-symmetric relax- 
ational models: (a) response propagator; (b) noise vertex; (c) anharmonic vertex. 



in Eq. H3fi|l as a four-point vertex. In the diagrams representing the different 
terms in the perturbation series, these vertices serve as hnks for the propa- 
gator hnes, with the fields 5" being encoded as the 'incoming', and the 
as the 'outgoing' components of the hnes. In Fourier space, translational in- 
variance in space and time imphes wavevector and frequency conservation at 
each vertex, see Fig. |21 belovif. An alternative, equivalent representation uses 
both the response and correlation propagators as independent elements, the 
latter depicted as undirected line, thereby disposing of the noise vertex, and 
retaining the nonlinearity in Fig. Q^c) as sole vertex. 

Following standard field theory procedures m^HHEE], one establishes 
that the perturbation series for the cumulants H41|l is given in terms of con- 
nected Feynman graphs only (for a detailed exposition of this and the following 
results, see Ref. ^Sj)- An additional helpful reduction in the number of dia- 
grams to be considered arises when one considers the vertex functions, which 
generalise the self-energy contributions S{q, lo) in the Dyson equation for the 
response propagator, G{q,uj)~^ — Dq"^ x{Qt'^)~^ = ^0(5,^)"^ — U{q,uj). 
To this end, we define the fields ^" = 6\nZ/5j°' and 'P" = S\nZ/Sj°', and 
introduce the new generating functional 

r[$, <P]^- lnZ[j, j] + fd'xfdtY, {r + r <^") , (43) 

•' a 

wherefrom the vertex functions are obtained via the functional derivatives 

_ N ^ N 



i j 



(44) 



Diagrammatically, these quantities turn out to be represented by the possible 
sets of one-particle ( IPI) irreducible Feynman graphs with N incoming and 
N outgoing 'amputated' legs; i.e., these diagrams do not split into allowed 
subgraphs by simply cutting any single propagator line. For example, for the 
two-point functions a direct calculation yields the relations 

r(i'i)(q,c.) = Dq'^xi-q,-^)-' = Go{-q,-u:)-' - S{-q,-u:) , (45) 
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where the second equation for r^^'^"> follows from the fluctuation-dissipation 
theorem (|13ll . Note that /^(°'^^(q, — vanishes because of causality. 




a 

Fig. 2. One-loop diagrams for (a) f'^'^' and (b) _r'^''^' in the time domain. 



The perturbation series can then be organised graphically as an expansion 
in successive orders with respect to the number of closed propagator loops. As 
an example, Fig. |21 depicts the one- loop contributions for the vertex functions 
r'(-'^d) and r^^'^^ in the time domain with all required labels. One may for- 
mulate general Feynman rules for the construction of the diagrams and their 
translation into mathematical expressions for the Ith. order contribution to 
the vertex function r^-^'-^"^: 

1. Draw Jill topologically diff'erent, connected one-particle irreducible graphs 
with N outgoing and N incoming lines connecting / relaxation vertices oc 
u. Do not allow closed response loops (since in the Ito calculus 0(0) = 0). 

2. Attach wavevectors qi, frequencies uji or times ti, and component indices 
ttj to all directed lines, obeying 'momentum (and energy)' conservation at 
each vertex. 

3. Each directed line corresponds to a response propagator Go(— g, — ti)) or 
Go('Z, ti—tj) in the frequency and time domain, respectively, the two-point 
vertex to the noise strength 2Dq'^, and the four-point relaxation vertex 
to —Dq°-u/6. Closed loops imply integrals over the internal wavevectors 
and frequencies or times, subject to causality constraints, as well as sums 
over the internal vector indices. Apply the residue theorem to evaluate 
frequency integrals. 

4. Multiply with —1 and the combinatorial factor counting all possible ways 
of connecting the propagators, I relaxation vertices, and k two-point ver- 
tices leading to topologically identical graphs, including a factor l//!fc! 
originating in the expansion of exp(— ^int[5', S]). 

For later use, we provide the explicit results for the two-point vertex func- 
tions to two-loop order. After some algebra, the three diagrams in Fig. Ogive 



r^'^'^\q,uj) =iuj + Dq'' 



, n+2 f 1 



■ u 



6 ; J^r + k^ Jk, [r + k'^ 
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Fig. 3. One-particle irreducible diagrams for r^^'^\q,Lu) to second order in 



n + 2 
18 



1 



1 



X 1 



kV + k^ J,,, r + k'^ r+{q-k- k' 

iuj 



iLu + A{k) + A{k') + A{q-k- k') 



(47) 



where we have separated out the dynamic part in the last line, and introduced 
the abbreviations A{q) = Dq" (r + q^) and = /d''/c/(27r)'* JJ^. For the 



(a) 





Fig. 4. (a) Two-loop diagram for r^'^'^\q,Lj); (b) one-loop graph for J^'^'' 



noise vertex, Fig. QJa) yields |E| 



1 



r + (q — fe — fe 



t\2 



l + Dq" 



Rc 



n + 2 



18 " .IkT + k-^jk, 
1 



r + k' 



iLo + A{k) -f A{k') + A{q-k- k') 



■ (48) 



notice that for model B, as a consequence of the conservation law for the 
order parameter and ensuing wavevector dependence of the nonlinear vertex, 
see Fig. ^c), to all orders in the perturbation expansion 



= 2 : r(i'i)(g = 0,c.)=ic., ^^r^^-^){q,u) 



= -2D 



(49) 



q=0 



At last, with the shorthand notation k = {q,uj), the analytical expression 
corresponding to the graph in Fig. 0{b) for the four-point vertex function at 
symmetrically chosen external wavevector labels reads 
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r^'-'\-^k/2;{k/2})=D{lq 



, n + 8 

„ ^ , „ 1 — u 

2^ J 6 

11/ icu 



Ik r + k'^ r + {q - k)'^ \ iuj + A{k) + A{q - k) 
1.4 Renormalisation 



(50) 



Consider a typical loop integral, say the correction in Eq. (|50ll to the four-point 
vertex function r^^'^^ at zero external frequency and momentum, whose 'bare' 
value, without any fluctuation contributions, is u. In dimensions < 4, one 
obtains, after introducing d-dimensional spherical coordinates and rendering 
the integrand dimensionless {x = \k\/^): 



■dx, (51) 



(27r)'' (r + fc2)2 2d-^nd/2r{d/2) Jo {I + x^) 

where we have inserted the surface area Sd = 2'k'^/'^ / r{d/2) of the d- 
dimensional unit sphere, with Euler's Gamma function, 7^(1 + x) — x r{x). 
Note that the integral on the right-hand side is finite. Thus, we see that the 
effective expansion parameter in perturbation theory is not just u, but the 
combinaton Ucff = Mr^'*"*^/^. Far away from Tc, it is small, and the pertur- 
bation expansion well-defined. However, Ueff ^ cxi as r ^ for d < 4: we 
are facing infrared (IR) divergences, induced by the strong critical fluctua- 
tions that render the loop corrections singular. A straightforward application 
of perturbation theory will therefore not provide meaningful results, and we 
must expect the fluctuation contributions to modify the critical power laws. 

Conversely, for dimensions d > 4, the integral in H51|l develops ultraviolet 
(UV) divergences as the upper integral boundary is sent to infinity (fc — |fc|). 

In lattice models, there is a finite wavevector cutoff, namely the Brillouin 
zone boundary, A ~ (27r/ao)'^ for a hypercubic lattice with lattice constant 
oo, whence physically these UV problems do not emerge. Yet we shall see that 
a formal treatment of these unphysical UV divergences will allow us to infer 
the correct power laws for the physical IR singularities associated with the 
critical point. The borderline dimension that separates the IR and UV singular 
regimes is referred to as upper critical dimension dc', here dc = 4. Note that 
at dc, UV and IR singularities are intimately connected and appear in the 
form of logarithmic divergences, see Eq. (|52|l . The situation is summarised 
in Table ^ where we have also stated that models with continuous order 
parameter symmetry, such as the Hamiltonian (0 with n > 2, do not allow 
long-range order in dimensions d < die = 2 {Mermin-Wagner-Hohenberg 
theorem |21l 1221 1^ ). Here, d\c is called the lower critical dimension; for the 
Ising model represented by Eq. Q with n = 1, of course die = 1. 
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Table 1. Mathematical and physical distinctions of the regimes d < dc, d — dc, and 
d > dc, for the 0(n)-symmetric models A and B (or static field theory). 



dimension 
interval 



perturbation 



model A / B or 
field theory 



critical 
behaviour 



2 < d < 4 



d = dc 



d > 4 



ill-defined no long-range 

u relevant order (n > 2) 

super-renormalisable nonclassical 
exponents 



d < die — 2 IR-singular 

UV-convergent 
IR-singular 

UV-convergent u relevant 
logarithmic IR-/ renormalisable 
UV-divergence u marginal 
IR-regular nonrenormalisable 
UV- divergent it irrelevant 



logarithmic 
corrections 
mean-field 
exponents 



The upper critical dimension can be obtained in a more direct manner 
through simple power counting. To this end, we introduce an arbitrary mo- 
mentum scale /i, i.e., define the scaling dimensions [x] = and [q] ~ fi. If in 
addition we choose [t] = or [us] = then the relaxation constant 

becomes dimensionless, [D] — jjP . For the deviation from the critical point, 
we obtain [r] — fj?, and the positive exponent indicates that this control pa- 
rameter constitutes a relevant coupling in the theory; as we shall see below, 
its renormahsed counterpart grows under subsequent RG transformations. For 
the nonlinear coupling, one finds [u] — /i'*^'*, so it is relevant for < 4: nonlin- 
ear thermal fluctuations will qualitatively affect the physical properties at the 
phase transition; but u becomes irrelevant for d > 4: one then expects mean- 
field (Gaussian) critical exponents. At the upper critical dimension dc — 4, 
the nonlinear coupling u is marginally relevant: this will induce logarithmic 
corrections to the mean- field scaling laws, see Tabled 

It is obviously not a simple task to treat the IR-singular perturbation 
expansion in a meaningful, well-defined manner, and thus allow nonanalytic 
modifications of the critical power laws (note that mean-field scaling is com- 
pletely determined by dimensional analysis or power counting). The key of 
the success of the RG approach is to focus on the very specific symmetry 
that emerges near critical points, namely scale invariance. There are several 
(largely equivalent) versions of the RG method; we shall here formulate and 
employ the field-theoretic variant 12 E El El El ■ order to pro- 
ceed, it is convenient to evaluate the loop integrals in momentum space by 
means of dimensional regularisation, whereby one assigns finite values even to 
UV-divergent expressions, namely the analytically continued values from the 
UV-finite range. For example, even for noninteger dimensions d and cr, we set 

f d'^k k^- _ ri<7 + d/2)r{s-a-d/2) 

J {2tt)'' {t + k^y 2dnd/2r{d/2)r{s) ' ^ ' 

The renormalisation program then consists of the following steps: 
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1. We aim to carefully keep track of formal, unphysical UV divergences. In 
dimensionally regularised integrals lf5S|l . these appear as poles in e = dc—d; 
their residues characterise the asymptotic UV behaviour of the field theory 
under consideration. 

2. Therefrom we may infer the (UV) scaling properties of the control pa- 
rameters of the model under a RG transformation, namely essentially a 
change of the momentum scale fi, while keeping the form of the action 
invariant. This will allow us to define suitable running couplings. 

3. We seek fixed points in parameter space where certain marginal couplings 
(u here) do not change anymore under RG transformations. This describes 
a scale-invariant regime for the model under consideration, where the 
UV and IR scaling properties become intimately linked. Studying the 
parameter flows near a stable RG fixed point then allows us to extract 
the asymptotic IR power laws. 

As a preliminary step, we need to take into account that the fluctuations 
will also shift the critical point downwards from the mean-field phase transi- 
tion temperature T°; i.e., we expect the transition to occur at Tc < T°. This 
fluctuation-induced Tc shift can be determined by demanding that the inverse 
static susceptibility vanish at Tc- x{q = 0,a; = 0)""'^ = r = r — Tc, where 
T ^ T -Tc and thus rc = Tc - T°. Using our previous results lHU and (jTzjl . 
we find to first order in u (and with finite cutoff yl), 

" + 2 /"^ 1 2x n + 2 uSdA''-^ ^. 2, ,,,, 



6 Jk rc + k"^ 6 (27r)'^(d-2) 

Notice that this quantity depends on microscopic details (the lattice structure 
enters the cutoff yl) and is thus not universal; moreover it diverges for d > 2 
(quadratically near dc = 4) as yl ^ oo. We next use r = r + Tc to write 
physical quantities as functions of the true distance t from the critical point, 
which technically amounts to an additive renormalisation; e.g., the dynamic 
response function becomes to one-loop order 



/ \— 1 '"^ 9 



n + 2 
1 u 



6 7fefc2(T + fc2 



0(^2). (55) 



The remaining loop integral is UV-singular in dimensions d > dc = 4:. 

We may now formally absorb the remaining UV divergences into renor- 
malised fields and parameters, a procedure called multiplicative renormalisa- 
tion. For the renormalised fields, we use the convention 

S% = 5" , 5^ = zl^^S" , (56) 

where we have exploited the 0{n) rotational symmetry in using identical 
renormalisation constants (Z factors) for each component. The renormalised 
cumulants with N order parameter fields 5" and N response fields S'" natu- 
rally involve the product Zg^"^ ^^^^ ' whence 
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4^-^) = Z-/^^ r(^-^) . (57) 

In a similar manner, we relate the 'bare' parameters of the theory via Z factors 
to their renormalised counterparts, which we furthermore render dimension- 
less through appropriate momentum scale factors, 

Dji^ ZdD , tr = ZrT , UB. ^ ZuuAd n'^^'^ , (58) 

where we have separated out the factor Ad — r{S — d/2)/2'^^^ n'^^^ for con- 
venience. In the minimal subtraction scheme, the Z factors contain only the 
UV-singular terms, which in dimensional regularisation appear as poles at 
e = 0, and their residues, evaluated at d = dc- 

These renormalisation constants are not all independent, however; since 
the equilibrium fluctuation-dissipation theorem H38() or (|13ll must hold in the 
renormalised theory as well, we infer that necessarily 

Zn=[Zs/Z~y^^ , (59) 
and consequently from Eq. H45|l 

Xb. = ZsX ■ (60) 

Moreover, for model B with conserved order parameter Eq. (|49|l implies that 
to all orders in the perturbation expansion 

a = 2 : Z~Zs = l , Zd = Zg . (61) 

For the following, it is crucial that the theory is renormalisable, i.e., a finite 
number of reparametrisations suffice to formally rid it of all UV divergences. 
Indeed, for the relaxational models A and B, and the static Ginzburg-Landau- 
Wilson Hamiltonian lO , all higher vertex function beyond the four-point func- 
tion are UV-convergent near dc, and there are only the three independent static 
renormalisation factors Zs, Zr^ and Z^, and in addition Zjj for nonconserved 
order parameter dynamics. As we shall see, these directly translate into the 
two independent static critical exponents and the unrelated dynamic scaling 
exponent z for model A; for model B with conserved order parameter, Eq. H61|) 
will yield a scaling relation between z and 77. 

In order to explicitly determine the renormalisation constants, we need 
to ensure that we stay away from the IR-singular regime. This is guaranteed 
by selecting as normalisation point either = 1 (i.e., ZrT = /i^) or q = /i. 
Inevitably therefore, the renormalised theory depends on the corresponding 
arbitrary momentum scale /i. Since there are no fluctuation contributions to 
order u to either dr'^^'^^q,0)/dq'^ or dr'^^'^\0,uj)/duj (at tr = 1), we find 
Zs = 1 and Zd = 1 within the one-loop approximation. Expressions 1)55(1 and 
then yield with the formula 
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Zr = 1 



n + 2 
~6 7 



1 



n + 8 w^j 

~6 r 



(62) 



To two-loop order, we may infer the field renormalisation Zg from the static 
susceptibility as the singular contributions to 9xfi(q, 0)/9q^|q=o, and Zd for 
model A, through a somewhat lengthy calculation [T^, from either r'^^'°''(0, 0) 
or I^jj^'^-' (0, cj), with the results 



Z. = \ 



n + 2ul n + 2 

— , a — : Zn = 1 

144 e ' 144 



6 In - - 1 



(63) 



1.5 Scaling laws and critical exponents 

We now wish to related the renormalised vertex functions at different inverse 
length scales /i. This is accomplished by simply recalling that the unrenor- 
malised vertex functions obviously do not depend on fi, 



d/i 



Z. 



N/2 yN/2 UN,N) 



(64) 



In the second step, the bare quantities have been replaced with their renor- 
malised counterparts. The innocuous statement H64|) then implies a very non- 
trivial partial differential equation for the renormalised vertex functions, the 
desired renormalisation group equation^ 



d Nj~ + N-fs 



d 



ItTR 



d 

dTf, 



■I3u 



d 



duf 



{N,N) 



{^i,D]i,tr,ur) = 



(65) 



Here we have defined Wilson's flow functions (the index '0' indicates that 
the derivatives with respect to fi are to be taken with fixed unrenormalised 
parameters) 



d 
d 
d 

1D = ^.-^ 



\nZ~ 



7s = M 



d_ 



InZs 



Htr/t) = -2 + fi — 
HDr/D) = J f 75 - 7 



In Zr , 



(66) 
(67) 
(68) 



where we have used the relation (|59|l : for model B, Eq. H61|l gives in addition 



7-D = 7s = -7^ 



(69) 



We have also introduced the RG beta function for the nonlinear coupling u, 
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R ^ 

on 



Ufl Ufl « - 4 + /I — 
V 



lnZ„ 



(70) 



Explicitly, Eqs. and H62I) yield to lowest nontrivial order, with e = 4 — d, 

n + 2 



IS 

a = : 7d 



72 

n + 2 



4 



72 >61n--l)ii^ + OK) , 



n + 2 

7r = -2 + — ^ + O(u^) 



+ 8 ^/ 2 N 



(71) 
(72) 
(73) 
(74) 



In the RG equation for the renormalised dynamic susceptibility, Eq. Ht)U|) 
tells us that the second term in Eq. H65|l is to be replaced with —75'. Its explicit 
dependence on the scale can be factored out via xnifJ-T DR^Tfi,uji,q,Lu) — 
1^'"^ XR {■tr,ur, q/ n,LO / Dr , see Eq. whence 



^ d d ^ d 

-2 - 7s + 7£) Dr — — + 7^ TR h Pu ^ — 

ODr dTR duR 



Xr{Dr,tr,ur) = . (75) 



This linear partial differential equation is readily solved by means of the 
method of characteristics, as is Eq. ll^ for the vertex functions. The idea 
is to find a curve parametrisation ^,{1) — ii£ in the space spanned by the 
parameters D, r, and u such that 

with initial values Dr, tr, and ur, respectively at £ = 1. The first- order or- 
dinary differential equations H76(l . with 7d(^) = 1d{u{C)^ etc. define running 

couplings that describe how the parameters of the theory change under scale 
transformations ji ^ fi£. The formal solutions for the D{£) and D(£) read 



D{£) ^ Dr exp 





df ■ 




' £' 



t{£) = TR exp 



^ d£' 



(77) 



For the function x{£) — xr yp{£), t{£), u{£) ) , we then obtain another ordinary 
differential equation, namely 



d£ 



(78) 



which is solved by 
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x{l)^xWe exp 



Collecting everything, we finally arrive at 



(79) 



XR{fJ.,DR,TR,UR,q,uj) = {p£) ^ cxp 



/ isii') 



2+a 



(80) 



The solution (|80|l of the RG equation (|75(l . along with the flow equations 
H76(l . (j77|l for the running couplings tell us how the dynamic susceptibility 
depends on the (momentum) scale /li^ at which we consider the theory. Sim- 
ilar relations can be obtained for arbitrary vertex functions by solving the 
associated RG equations (|65|) ^13^ . The point here is that the right-hand side 
of Eq. (|5n|) may be evaluated outside the IR-singular regime, by fixing one 
of its arguments at a finite value, say \q\/^(. = 1. The function xr is regu- 
lar, and can be calculated by means of perturbation theory. A scale- invariant 
regime is characterised by the renormalised nonlinear coupling ur becoming 
independent of the scale or m(^) — > u* = const. For an RG fixed point to 
be infrared- stable^ we thus require 



/?„(u*)=0, /?:(«*) >0 



(81) 



since Eq. (|7()|l then implies that u{l — > 0) ^ m*. Taking the hmit i ^ Q 
thus provides the desired mapping of physical observables such as H80(l onto 
the critical region. In the vicinity of an IR-stable RG fixed point, Eq. H77(l 
yields the power laws D{(.) « DrP'o^ where 7^ = 7d(^ ^ 0) = 7d(u*), etc. 
Consequently, Eq. l|Sn|) reduces to 



Xr{tr, q, u;) « r^^'s xr [tr l< , 



(82) 



and upon matching I = |q|//i we recover the dynamic scaling law \12\ with 
the critical exponents 

V^-l*s, v^-\hl, ^ = 2 + a + 7|5. (83) 
To one-loop order, we obtain from the RG beta function (|74|l 

6e 



iff 



n-l-8 



+ 0(6^) 



(84) 



Here we have indicated that our perturbative expansion for small u has ef- 
fectively turned into a dimensional expansion in e = dc — d. In dimensions 
d < 4, the Heisenberg fixed point u*^ is IR-stable, since f3'^{u*fj) = e > 0. With 
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Eqs. H71|l and O, the identifications H83|) then give us explicit results for the 
static scaling exponents, as mere functions of dimension = 4 — e and the 
number of order parameter components n, 

For model A with nonconserved order parameter, the two-loop result H72|l 
yields the independent dynamic critical exponent 

4 

a = : z^2 + c'n , c = 6 In - - 1 + 0(e) ; (86) 

o 

for model B with conserved order parameter, instead — Js ~ whence 
we arrive at the exact scaling relation 

a = 2 : z = 4- 7?. (87) 

In dimensions d > dc = A, the Gaussian fixed point Mq = is stable 
(/3(j(0) — — e > 0). Therefore all anomalous dimensions disappear, i.e., 7^ = 
= 7^ and 7* = —2, and we are left with the mean- field critical exponents 
?7o = 0, I'o = 1/2, and zq = 2 + a. Precisely at the upper critical dimension 
dc = 4, the RG flow equation for the nonlinear coupling becomes 



which is solved by 



^^ = ^W + 0(iI(f)3), (33) 



= ^ . (89) 



In four dimensions, u{l) 0, but only logarithmically slowly, which causes 
logarithmic corrections to the mean- field critical power laws. For exam- 
ple, upon inserting Eq. into the flow equation (f7?)|l . one finds t{£) ~ 
Ti? £-2(ln |^|)-("+2)/(«+8). ^ith T{e = ^-1) = 0(1), iterative inversion yields 

e(r^)^r^^/^(h.r«)("+^)/^("+«) . (90) 

This concludes our derivation of asymptotic scaling laws for the critical 
dynamics of the purely relaxational models A and B, and the explicit com- 
putation of the scaling exponents in powers of e = dc — d. In the following 
sections, I will briefly sketch how the response functional formalism and the 
dynamic renormalisation group can be employed to study the critical dynam- 
ics of systems with reversible mode-coupling terms, the 'ageing' behaviour 
induced by quenching from random initial conditions to the critical point, the 
effects of violating the detailed balance constraints on universal dynamic crit- 
ical properties, and the generically scale-invariant features of nonequihbrium 
systems such as driven diffusive Ising lattice gases. 
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1.6 Critical dynamics with reversible mode-couplings 

In the previous chapters, we have assumed purely relaxational dynamics for 
the order parameter, see Eq. H16|l . In general, however, there are also re- 
versible contributions to the systematic force terms F"' that enter its Langevin 
equation [7l l24j . Consider the Hamiltonian dynamics of microscopic variables, 

say, local spin densities, at T = 0: dtS'^{x,t) — |_ff[S'm], 5*^(2;, Here, 

the Poisson brackets {A, B} constitute the classical analog of the quantum- 
mechanical commutator B] (correspondence principle). Upon coarse- 
graining, the microscopic variables become the mesoscopic hydrodynamic 
fields 5" . Since the set of slow modes should provide a complete description 
of the critical dynamics, we may formally expand 



(91) 



{H[^],^"(.)}=|dvx:^Q^"(-',-), 

with the mutual Poisson brackets of the hydrodynamic variables 

Q"^ix, x') = {s"{x),S^ix')} = -Q^'^ix', x) . (92) 

By inspection of the associated Fokker-Planck equation, one may then 
establish an additional equilibrium condition in order for the time-dependent 
probability distribution to reach the canonical limit ©: 7^[5, t] VcqiS] as 
t —^ oo provided the probability current is divergence-free in the space spanned 
by the stochastic fields S°'{x): 

/d^-E^(^."v[^]e--[-l/^--)^0. (93) 

It turns out that this equilibrium condition is often more crucial than the 
Einstein relation (|17|l . In order to satisfy Eq. H93|) at T ^ 0, we must sup- 
plement Eq. (|91|l by a finite-temperature correction, whereupon the reversible 
mode- coupling contributions to the systematic forces become 



and the complete coupled set of stochastic differential equations reads 



(94) 



= F,ZAS]{x,t) - i^"(iV)-^|^ 4- C{x,t) , (95) 

where as before the D" denote the relaxation coefficients, and a" = or 2 
respectively for nonconserved and conserved modes. 
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isotropic ferromagnets, H[{Sj}] — — ^ X]ffc=i '^jk '^j ' '^k, where the spin op- 



As an instructive example, let us consider the Heisenberg model for 

1 ^AT 

erators satisfy the usual commutation relations 

The corresponding Poisson brackets for the magnetisation density read 



''^{x,x') = -gY^e'^^-'S''ix)S{x-x') , (96) 



where the purely dynamical coupling g incorporates various factors that 
emerge upon coarse-graining and taking the continuum limit. The second 
contribution in Eq. (|94(l vanishes, since it reduces to a contraction of the an- 
tisymmetric tensor e"'^'*' with the Kronecker symbol 5^^'* , whence we arrive 
at the Langevin equations governing the critical dynamics of the three order 
parameter components for isotropic ferromagnets |25| 

ds{x,t) ^, , 5n[s] ^ . 5n[s] ^, , 

with (C(a?,t)) — 0. Since ii\{S jW'Yl,k^k — 0' ^he total magnetisation is 
conserved, whence the noise correlators should be taken as 

(C"(a;,i)C'^(a;',t')) = ~2Dk^TV'^8{x- x')8{t~t')&°'^ . (98) 

The vector product term in Eq. (|97|l describes the spin precession in the local 
effective magnetic field (57i[5']/(55, which includes a contribution induced by 
the exchange interaction. 

The Langevin equation (|97(l and H98() with the Hamiltonian ((Jj) for n = 3 
define the so-called model J 7 . In addition to the model B response functional 
H35|) and H36|l with a — 2 (setting k-oT = 1 again), the reversible force in 
Eq. (|97|l leads to an additional contribution to the action 



j d'^x J<itYl ^"^^S^S^ {V^S^ + h'') , (99) 



which gives rise to an additional mode-coupling vertex, as depicted in Fig.|S{a). 
Power counting yields the scaling dimension [g] = /i^"''/^ for the associated 
coupling strength, whence we expect a dynamical upper critical dimension 
d'^ — 6. However, since we are investigating a system in thermal equilibrium, 
we can treat its thermodynamics and static properties separately from its 
dynamics. Obviously therefore, the static critical exponents must still be given 
(to lowest nontrivial order and for d < dc = 4) by Eqs. H85|l for the three- 
component Heisenberg model with 0(3) rotational symmetry. Therefore our 
sole task is to find the dynamic critical exponent z. 

Remarkably, z is entirely fixed by the symmetries of the problem and can 
be determined exactly. To this end, we exploit the fact that the 5" are the 
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(a) (b) (c) 



1 




Fig. 5. (a) Mode-coupling three-point vertex for model J. One-loop Feynman dia- 
grams for the propagator (b) and noise vertex (c) renormalisations in model J. The 
same graphs (b), (c) apply for driven diffusive systems (Sec. 11.91 . 



generators of the rotation group; indeed, it follows from Eq. (|99|l that applying 
a time-dependent external field h^it) induces a contribution 

S"{x,t)) =g fdt'"S2e''>^''(s^{x,t')) h-^{t) (100) 

to the average magnetisation. As a consequence, we obtain for the nonlinear 
susceptibility i?"''^')' = 6^{S°')/Sh'^ <5/i''U=o, 

j d'^x' R°'-'^^{x,t;x~ x',t-t') = g e"'^^ x'^Hx^t) 0{t) 0{t - t') . (101) 

An analogous expression must hold after renormalisation as well. If we define 
the dimensionless renormalised mode-coupling according to 

gl^Zgg^Ba/-' , fR^gl/Dl, (102) 

where Bd = r{4 - d/2)/2'^ dn'^/'^ , Eq. (fTnT|l implies the identity 

Zg = Zs . (103) 

For the RG beta function associated with the effective coupling entering the 
loop corrections, we thus infer 



ft ^ 



fB = fRid-6 + ls-2jD) ■ (104) 





Consequently, at any nontrivial IR-stable RG fixed point < /* < cx3, we 
have the exact scaling relation, valid to all orders in perturbation theory, 

d<6 : z = 4 + 7|,=4-f ^ = (105) 

Since the resulting value for the dynamic exponent, z « 5/2 in three dimen- 
sions, is markedly smaller than the model B mean-field zq = 4, we conclude 
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that the reversible spin precession kinetics speeds up the order parameter 
dynamics considerably PSIITIIU]. 

An explicit one-loop calculation, either for the propagator self-energy 
(q, w), depicted in Fig. [S{b), or the noise vertex r(2.o)(qr,a;), shown m 
Fig.lEc), yields Ii[T3| 

7D = -fB. + 0{ulJl,) , (106) 

which along with 7s = -I- 0(w|j, /^) confirms that there exists a nontrivial 
mode-coupling RG fixed point 

where e = 6 — d, which is IR-stable for d < 6. As ?7 = for d > 4, we 
indeed recover the mean- field dynamic exponent zq = 4 in d > 6 dimensions. 
With the leading singularity thus isolated, the regular scaling functions can 
be computed numerically to high accuracy within a self-consistent one-loop 
approximation that also goes under the name mode- coupling theory. Details of 
this procedure, an alternative derivation, and many results of mode-coupling 
theory as applied to the critical dynamics of magnets and comparisons with 
experimental data can be found in Ref . 26 . 

Typically, reversible force terms of the form 194(1 involve dynamical cou- 
plings of the order parameter to other conserved, slow variables. In addition, 
there may also be static couplings to conserved fields in the Hamiltonian. 
These various possibilities give rise to a range of different dynamic universal- 
ity classes for near-equilibrium critical dynamics 7 . We shall not pursue these 
further here (for a partial account within the field-theoretic RG approach, see 
Ref. jl8|l. but instead proceed and now consider nonequihbrium effects. 



1.7 Critical relaxation, initial slip, and ageing 

We begin with a brief discussion of the coarsening dynamics of systems de- 
scribed by model A / B kinetics that are rapidly quenched from a disordered 
state at T > Tc to the critical point T k E71EHI- The situation may 
be modeled as a relaxation from Gaussian random initial conditions, i.e., the 
probability distribution for the order parameter at t = can be taken as 

P[5',t = 0] cxe^^«[^l =exp ^-^y"d'*x^[S'"(a;,0)-a"(a;)]^^ , (108) 

where the functions a"(£c) specify the most likely initial configurations. Power 
counting for the parameter A gives [A] = whence it is a relevant pertur- 
bation that will flow to Z\ ^ cx) under the RG. Asymptotically, therefore, the 
system will be governed by sharp Dirichlet boundary conditions. Whereas the 
response propagators remains a causal function of the time difference between 
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applied perturbation and effect, Go{q,t - t') ^ 0{t - f) 6-'^''?° (*"*'), 
see Eq. (|23ll . time translation invariance is broken by the initial state in the 
Dirichlet correlator of the Gaussian model, 

Cj,{q; t, t') = { 0-^"° - e-^-J" . (109) 

r + V / 

Away from criticality, i.e., for r > and q ^ 0, temporal correlations decay 
exponentially fast, and the system quickly approaches the stationary equilib- 
rium state. However, as T ^ Tc, the equilibration time diverges according to 
tc ~ l''"!"^'' ~^ and the system never reaches thermal equilibrium. Two- 
time correlation functions will then depend on both times separately, in a 
specific manner to be addressed below, a phenomenon termed critical 'ageing ' 
(for more details, see Refs. [211 ED])- 

The field-theoretic treatment of the model A / B dynamical action (|35|l . 
I)3H|I with the initial term H108(l follows the theory of boundary critical phenom- 
ena |31| . However, it turns out that additional singularities on the temporal 
'surface' at i + t' = appear only for model A, and can be incorporated into 
a single new renormalisation factor; to one-loop order, one finds [271 128| 

a = : §g(a;,0) = (ZoZ~)i/2 5"(a;,0) , Zq = 1 - ^ ^ . (110) 

This in turn leads to a single independent critical exponent associated with 
the initial time relaxation, the initial slip exponent, which becomes for the 
purely relaxational models A and B with nonconserved and conserved order 
parameter: 

a.O: ..|.^. + 0(.^), a = 2:e^0. (Ill) 

In order to obtain the short-time scaling laws for the dynamic response and 
correlation functions in the ageing limit t' /t — > 0, one requires additional 
information that can be garnered from the short-distance operator product 
expansion for the fields, 

t^Q : S{x,t)^d{t)So{x) , S{x,t) = a{t)Sa{x) . (112) 

Subsequent analysis then yields eventually j57| 1^ 

x(q;t,t'-^0) 

C(q;t,t'^0) 

and for the time dependence of the mean order parameter 



= 191-'+" 



Co{qi,\q\'Dt) , 



(113) 
(114) 
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(Sit)) = Soto's [Sot'>'+^/''') , (115) 

a = : e' = 6- ^^^^ , = 2:6*' = 61 = 0. (116) 

z 

One may also compute the universal fluctuation-dissipation ratios in this 
nonequilibrium ageing regime [291 13U| . It emerges, though, that these depend 
on the quantity under investigation, which prohibits a unique definition of an 
effective nonequilibrium temperature for critical ageing. The method sketched 
above can be extended to models with reversible mode-couplings '32\ For 
model J capturing the critical dynamics of isotropic ferromagnets, one finds 

n_ + V _ 6-d~7J _ ^^^^^ 



z 

in systems where a nonconserved order parameter is dynamically coupled to 
other conserved modes, the initial slip exponent 9 is actually not a universal 
number, but depends on the width of the initial distribution i32| . 



1.8 Nonequilibrium relaxational critical dynamics 

Next we address the question What happens if the detailed balance 
conditions (|17|l and H93|l are violated? To start, we change the noise strength 
D D in the purely relaxational models A and B, which (in our units) 
violates the Einstein relation (|17|l . However, this modification can obviously be 
absorbed into a rescaled effective temperature^ k^T k^T' = D/D. Formally 
this is established by means of the dynamical action 134() . which now reads 



A[S,S] = J^'^xJ dtJ2S" 



dt S"' + D {iVy (r- V^) 5" 



D (iVf S'^+D^ {iVf S^S^ 



(118) 



Upon simple rescaling S*" ^ 5 " = S"" \Jd/D, S"' -> S'" = \Jd/D, the 
response functional H118(l recovers its equilibrium form, albeit with modified 
nonlinear coupling u — > u = uD/D. However, the universal asymptotic prop- 
erties of these models are governed by the Heisenberg fixed point (|84|l . and the 
specific value of the (renormalised) coupling, which only serves as the initial 
condition for the RG flow, does not matter. In fact, the relaxational dynamics 
of the kinetic Ising model with Glauber dynamics (model A with n = 1) is 
known to be quite stable against nonequilibrium perturbations |34[ , even 
if these break the Ising Z2 symmetry j^S]- For model J the above rescaling 
modifies in a similar manner merely the mode-coupling strength in Eq. (jSHJ, 



namely g g = g \/ D / D . Again, since the dynamic critical behaviour 
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is governed by the universal fixed point H1U7|) . thermal equilibrium becomes 
effectively restored at criticality. More generally, it has been established that 
isotropic detailed balance violations do not affect the universal properties in 
other models for critical dynamics that contain additional conserved variables 
either: the equilibrium RG fixed points tend to be asymptotically stable 

In systems with conserved order parameter, however, we may in addition 
introduce spatially anisotropic violations of Einstein's relation; for example, in 
model B one can allow for anisotropic relaxation —D —D± — V| , 

with different rates in two spatial subsectors and concomitantly anisotropic 
noise correlations — D —D± — V|. We have thus produced a 

truly nonequilibrium situation provided D±/D± ^ D^^/D^^, which we may 
interpret as having effectively coupled the longitudinal and transverse spatial 
sectors to heat baths with different temperatures T± < Ty, say |38| . 

Evaluating the fluctuation-induced shift of the transition temperature, 
see Eq. (j54|l one finds not surprisingly that the transverse sector softens 
first, while the longitudinal sector remains noncritical. This suggests that we 
can neglect the nonlinear longitudinal fluctuations as well as the Vjj term 
in the propagator. These features are indeed encoded in the corresponding 
anisotropic scaling: [q±] = /i, [gy] = /it^, [uj] = whence [D±] ~ [D±] = /i", 

and [D||] = [D||] ~ ji^^ become irrelevant. Upon renaming D = D±^ and 
c = r||_Dy/D^, this ultimately leads to the randomly driven or two-temperature 
model B [351 1301 as the effective theory describing the phase transition: 

dS"ix,t) _^ 
dt 



Vl{r-Vl)+cV\\S^[x,t) 

n 

6" 



+ ^VlS-{x,t)Y,[S^{x,t)r + C{x,t) , (119) 



13 

with the noise correlations 

{Cix, t) C^{x', t')) ^ ^2D Vl Six - x') Sit - t') 6"^ . (120) 

Quite remarkably, the Langevin equation (|119|l can be derived as an equilib- 
rium diffusive relaxational kinetics 

dS'^ix,t) _ sn^s] 



dt ^ 6S°'{x,t) 

from an effective long-range Hamiltonian 



D^i .or^: (121) 



- / 2. 1^ ( 



Jd''xj2[s'^ix)r[sP{x)r . (122) 
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Power counting gives [u\ — //^^''ii^'': the spatial anisotropy suppresses longi- 
tudinal fluctuations and lower the upper critical dimension to dc = 4 — d|| . 
The anisotropic correlations encoded in Eq. H122|l also reduce the lower critical 
dimension and affect the nature of the ordered phase j39l HT] . 

The scaling law for, e.g., the dynamic response function takes the form 

x(r,.,„„.„)H,x|-»'x(|^,|^,j^) , (123) 

where we have introduced a new anisotropy exponent A. Since the nonlinear 
coupling u only affects the transverse sector, we find to all orders in the 
perturbation expansion: 

r^'''\q^ = 0,q\\,u;) = iiu + D cql , (124) 

and consequently obtain the Z factor identity 

Zc = = Zg' , (125) 

which at any IR-stable RG fixed point implies the exact scaling relations 

z = 4-^, Zi=l-^ = l-| = £-l, (126) 

whereas the scaling exponents for the longitudinal sector read 

^|| = Y^ = 2> ^^11 =;^(l + /A) = ^(4-7?) . (127) 

As for the equilibrium model B, the only independent critical exponents to 
be determined are 77 and u. To one-loop order, only the combinatorics of the 
Feynman diagrams (see Fig. ^ enters their explicit values, whence one finds 
for d < dc = 4 — d|| formally identical results as for the usual Ginzburg- 
Landau- Wilson Hamiltonian 

V = + O{e'), ^=2-^e + 0i^), (128) 
I' n + 8 

albeit with different e — 4 — d — To two- loop order, however, the anisotropy 
manifestly affects the evaluation of the loop contributions, and the value for 
77 deviates from the expression in Eq. H85|l 40 . 

Interestingly, an analogously constructed nonequilibrium two-temperature 
model J with reversible mode-coupling vertex cannot be cast into a form that 
is equivalent to an equilibrium system, for owing to the emerging anisotropy, 
the condition cannot be satisfied. A one-loop RG analysis yields a run- 
away flow, and no stable RG fixed point is found |38| . Similar behaviour ensues 
in other anisotropic nonequilibrium variants of critical dynamics models with 
conserved order parameter; the precise interpretation of the apparent insta- 
bility is as yet unclear [6'6\ . 
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1.9 Driven diffusive systems 

Finally, we wish to consider Langevin representations of genuinely nonequi- 
librium systems, namely driven diffusive lattice gases (for a comprehensive 
overview, see Ref. j42j). First we address the coarse-grained continuum version 
of the asymmetric exclusion process, i.e., hard-core repulsive particles that hop 
preferentially in one direction. We describe this system in terms of a conserved 
particle density, whose fluctuations we denote with ^(a;, t), such that (S) = 0, 
obeying a continuity equation dt S{x, t) -I- V • J(a;, t) = 0. We assume the sys- 
tem to be driven along the '||' direction; in the transverse sector (of dimension 
d± = d — 1) we thus just have a noisy diffusion current J± — —D Vj^S* -I- r], 
whereas there is a nonlinear term, stemming from the hard-core interactions, 
in the current along the direction of the external drive, with Jq || — const.: J|| — 
Jq II — £) c VyS* — i Dg 5*^-1-7711 . For the stochastic currents, we assume Gaussian 
white noise {r]i) = = (?7||) and {r]i{x, t) r]j{x', t')) = 2D S{x — x') 5{t — t') Sij, 
(r^ii {x, t) 7^11 (a;', t')) — 2Dc5{x — x') S{t — t'). Notice that since we are not in 
thermal equilibrium, Einstein's relation need not be fulfilled. We can however 
always rescale the field to satisfy it in the transverse sector; the ratio w = c/c 
then measures the deviation from equilibrium. These considerations yield the 
generic Langevin equation for the density fluctuations in driven diffusive sys- 
tems (DDS) mim 

= D (Vi + cVl) Six, t) + ^ V||5(a3, t)' + Cix, t) , (129) 

with conserved noise C = — V^-r/ — V||?7||, where (C) = and 

{Cix, t) C(a;', t')) = ~2D {y\ + ^ Vjj) Six - x') 5{t - t') . (130) 

Notice that the drive term cx g breaks both the system's spatial reflection 
symmetry and the Ising Z2 symmetry S — > —S. 

The corresponding Janssen-De Dominicis response functional (|34() reads 



dt 

+D (vi -Hcvf) 5 - ^ V|| 52 



(131) 



It describes a 'massless ' theory, hence we expect the system to be generically 
scale-invariant, without the need to tune it to a special point in parameter 
space. The nonlinear drive term will induce anomalous scaling in the drive 
direction, different from ordinary diffusive behaviour. In the transverse sector, 
however, we have to all orders in the perturbation expansion simply 

r^^^^\q^,q\\^Q,u)^iu + Dql, r'^^'^Hq^, q^^ = 0, to) = -2D ql , (132) 
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since the nonlinear three-point vertex, which is of the form depicted in 
Fig.[SJa), is proportional to iq^^. Consequently, 

Zs ^ Zd^I , (133) 

which immediately implies 

?7 = , z = 2 . (134) 

Moreover, the nonlinear coupling g itself does not renormalise either as a 
consequence of Galilean invariance. Namely, the Langevin equation (|129|l and 
the action H131|) are left invariant under Galilean transformations 

S'{x'j^,x\^,t') = S'(a;_L,2:|| - Dgvt,t) - v ; (135) 

thus, the boost velocity v must scale as the field S under renormalisation, and 
since the product Dgv must be invariant under the RG, this leaves us with 

Zg = ^5 ' = 1 • (136) 

The effective nonlinear coupling governing the perturbation expansion in 
terms of loop diagrams turns out to be g^/c^/^; if we define its renormalised 
counterpart as 

VR = Z^J^vCd^i''-\ (137) 

with the convenient choice Cd = r{2 — d/2) ^ we see that the asso- 
ciated RG beta function becomes 

l3, = VR(^d~2~^^)j . (138) 

At any nontrivial RG fixed point < t;* < oo, therefore 7* = |(d — 2). 
We thus infer that below the upper critical dimension dc — 2 for DDS, the 
longitudinal scaling exponents are fixed by the system's symmetry [431 144j , 

^ = -^ = ^, .„^^ = ^. (139) 

An explicit one-loop calculation for the two-point vertex functions, see 
Fig-EIb) and (c), yields 

7c = (3 + ^r) , 7£ = (3^,^' + 2 + 3wr) , (140) 
f3w = WR (7£ - 7c) = [wR - 1) [wR - 3) . (141) 

This establishes that in fact the fixed point w* = 1 is IR-stable (provided < 
V* < 00), which means that asymptotically the Einstein relation is satisfied 
in the longitudinal sector as well |43| . 
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In this context, it is instructive to make an intriguing connection with the 
noisy Burgers equation j45| . describing simphfied fluid dynamics in terms of 
a velocity field u{x, t): 



+ :^ V [uix, tr] = D vM^. t) + c{x, t) , 



(142) 



(G)=0, (Qix,t)Qix',t')) =-2DV,V,S{x-x')6{t-t') . (143) 

For Dg ~ 1, the nonlinearity is just the usual fluid advection term. In one di- 
mension, the Burgers equation (|142|l becomes identical with the DDS Langevin 
equation (|129l) . so we immediately infer its anomalous dynamic critical ex- 
ponent Z|| — "ijl. At least in one dimension therefore, it should represent 
an equilibrium system which asymptotically approaches the canonical distri- 
bution ©, where the Hamiltonian is simply the fluid's kinetic energy (and 
we have set fceT — 1). So let us check the equilibrium condition H93|l with 
Vcq[u] cx exp [-i / uixf d'^x] : 



■ |_Vit(a;, t) \ e ^ J ^ ' 



6u{x, t) 

= y" [2V • u{x, t) - u{x, t) ■ Vu{x, tf] d-^x e"* /«("'^*)' d''.' 

With appropriate boundary conditions, the first term here vanishes, but the 
second one does so only in d = 1: — j u (du^ / dx) dx — jv?(duldx)dx = 
^ J {du^ /dx) dx = Q. Driven diffusive systems in one dimension arc therefore 
subject to a 'hidden' fluctuation-dissipation theorem. 

To conclude this part on Langevin dynamics, let us briefly consider the 
driven model B or critical DDS |42| , which corresponds to a driven Ising lattice 
gas near its critical point. Here, a conserved scalar field 5* undergoes a second- 
order phase transition, but similar to the randomly driven case, again only the 
transverse sector is critical. Upon adding the DDS drive term from Eq. H129|l 
to the Langevin equation 1119(1 . we obtain 



dS{x,t) 
dt 



= D 



Vi (r 

Dg 



vi) +cvf 



S{x,t) + —VlS{x,tf 

D 



V\\S{x,tY + Q{x,t) 



(144) 



with the (scalar) noise specified in Eq. (|120|l . The response functional thus 
becomes 



A[S, S] 



d^x / dt S 



dS 



Vi (r 



S 



+D 



(145) 
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Power counting gives [g'^] = fJ.^^'^, so the upper critical dimension here is 
dc — 5, and [u\ = fi^~'^. The nonhnearity cx u is thus irrelevant and can 
be omitted if we wish to determine the asymptotic universal scaling laws; 
but recall that it is responsible for the phase transition in the system. The 
remaining vertex is then proportional to \q\\ , whence Eqs. I)132|l and H133() hold 
for critical DDS as well, and the transverse critical exponents are just those 
of the Gaussian model B, 

rj^O, z = 4. (146) 

In addition, Galilean invariance with respect to Eq. H135() and therefore 
Eq. (|136|l hold as before. With the renormalised nonlinear drive strength de- 
fined similarly to Eq. (|137|l . but a different geometric constant and the scale 
factor /x''"^, the associated RG beta function reads 



l3v^vR\d-b- - -icj , (147) 

which again allows us to determine the longitudinal scaling exponents to all 
orders in perturbation theory, for c? < dc = 5, 

7r 8-d 4 12 , 

It is worthwhile mentioning a few marked differences to the two-temperature 
model B discussed in Sec. ll.8l In DDS, there are obviously nonzero three-point 
correlations, and in the driven critical model B the upper critical dimension is 
dc = 5 as opposed to dc = 4 — d|| for the randomly driven version. Notice also 
that the latter is characterised by nontrivial static critical exponents, but the 
kinetics is purely diffusive along the drive direction, Z|| = 2. Conversely for the 
driven model B, only the longitudinal scaling exponents are non-Gaussian. 



2 Reaction— Diffusion Systems 

We now turn our attention to stochastic interacting particle systems, whose 
microscopic dynamics is defined through a (classical) master equation. Below, 
we shall see how the latter can be mapped onto a stochastic quasi- Hamiltonian 
in a second-quantised bosonic operator representation Taking the 

continuum on the basis of coherent-state path integrals then yields a field 
theory action that may be analysed by the very same RG methods as described 
before in Sees. 11.3111.51 ffor more details, see the recent overview [T^). 

2.1 Chemical reactions and population dynamics 

Our goal is to study systems of 'particles' A,B, . . . that propagate through 
hopping to nearest neighbors on a d-dimensional lattice, or via diffusion in 
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the continuum. Upon encounter, or spontaneously, with given stochastic rates, 
these particles may undergo species changes, annihilate, or produce offspring. 
At large densities, the characteristic time scales of the kinetics will be governed 
by the reaction rates, and the system is said to be reaction-limited. In contrast, 
at low densities, any reactions that require at least two particles to be in 
proximity will be diffusion-limited: the basic time scale will be set by the 
hopping rate or diffusion coefhcient. 

As a first approximation to the dynamics of such 'chemical' reactions, let us 
assume homogeneous mixing of each species. We may then hope to be able to 
capture the kinetics in terms of rate equations for each particle concentration 
or mean density. Note that such a description neglects any spatial fluctuations 
and correlations in the system, and is therefore in character a mean-field 
approximation. As a first illustration consider the annihilation oi k — I > 
particles of species A in the irreversible fcth-order reaction fc A — > I A, with rate 
A. The corresponding rate equation employs a factorisation of the probability 
of encountering k particles at the same point to simply the fcth power of the 
concentration a(t), 

a{t) = -(fc- /) Aa(i)'= . (149) 
This ordinary differential equation is readily solved, with the result 

k = l : a{t) = a(0)e-^* , (150) 

k>2 : a{t) = [a(0)i-^- + (k - l){k -l)\t\ '^'^''-^^ . (151) 

For simple 'radioactive' decay {k — 1), we of course obtain an exponential 
time dependence, as appropriate for statistically independent events. For pair 
(fc = 2) and higher-order {k > 3) processes, however, we find algebraic long- 
time behaviour, a{t) {Xt)~^^^''~^\ with an amplitude that becomes inde- 
pendent of the initial density a(0). The absence of a characteristic time scale 
hints at cooperative effects, and we have to ask if and under which circum- 
stances correlations might qualitatively affect the asymptotic long-time power 
laws. For according to Smoluchowski theory |12| . we would expect the anni- 
hilation reactions to produce depletion zones in sufficiently low dimensions 
d < dc, which would in turn induce a considerable slowing down of the den- 
sity decay, see Sec. 12.31 For two-species pair annihilation A-\- B ^ (i) (without 
mixing), another complication emerges, namely particle species segregation in 
dimensions for d < ds', the regions dominated by either species become largely 
inert, and the annihilation reactions are confined to rather sharp fronts |12| . 

Competition between particle decay and production processes, e.g., in the 
reactions A ^ $ (with rate k), A ^ A -\- A (with forward and back rates a 
and A, respectively), leads to even richer scenarios, as can already be inferred 
from the associated rate equation 

d{t) = {a - K)a{t) - Xa{tf . (152) 

For a < K, clearly a{t) ~ Q-i'-^-'^)^ ^ q as i ^ 00. The system eventually 
enters an inactive state, which even in the fully stochastic model is absorbing, 
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since once there is no particle left, no process whatsoever can drive the system 
out of the empty state again. On the other hand, for cr > k, we encounter an 
active state with a{t) — > = {a — k)/X exponentially, with rate ^ a — k. We 
have thus identified a nonequihbrium continuous phase transition at CTc = k. 
Indeed, as in equilibrium critical phenomena, the critical point is governed by 
characteristic power laws; for example, the asymptotic particle density aca ~ 
{a—GcY, and the critical density decay a{t) ~ (At)~" with /3o = 1 = ao in the 
mean-field approximation. The following natural questions then arise: What 
are the critical exponents once statistical fluctuations are properly included 
in the analysis? Can we, as in equilibrium systems, identify and characterise 
certain universality classes, and which microscopic or overall, global features 
determine them and their critical dimension? 

Already the previous set of reactions may also be viewed as a (crude) model 
for the population dynamics of a single species. In the same language, we may 
also formulate a stochastic version of the classic Lotka-Volterra predator-prey 
competition model j46| : if by themselves, the 'predators' A die out according 
to j4 — !■ 0, with rate k, whereas the prey reproduce B ^ B + B with rate a, 
and thus proliferate with a Malthusian population explosion. The predators 
are kept alive and the prey under control through predation, here modeled as 
the reaction A + B ^ A + A: with rate A, a prey is 'eaten' by a predator, who 
simultaneously produces an offspring. The coupled kinetic rate equations for 
this system read 

a{t) = Xa{t)b{t)- Ka{t) , h{t) ^ a b{t) - X a{t) b{t) . (153) 

It is straightforward to show that the quantity if (f) = X[a{t)+b{t)] — a\na{t)~ 
Khib{t) is a constant of motion for this coupled system of differential equa- 
tions, i.e., K{t) = 0. As a consequence, the system is governed by regular 
population oscillations, whose frequency and amplitude are fully determined 
by the initial conditions. Clearly, this is not a very realistic feature (albeit 
mathematically appealing), and moreover Eqs. (| 15 3(1 are known to be quite 
unstable with respect to model modifications 46 . Indeed, if one includes 
spatial degrees of freedom and takes account of the full stochasticity of the 
processes involved, the system's behaviour turns out to be much richer '47': 
In the species coexistence phase, one encounters for sufficiently large values of 
the predation rate an incessant sequence of 'pursuit and evasion' waves that 
form quite complex dynamical patterns. In finite systems, these induce erratic 
population oscillations whose features are however independent of the ini- 
tial configuration, but whose amplitude vanishes in the thermodynamic limit. 
Moreover, if locally the prey 'carrying capacity' is limited (corresponding to 
restricting the maximum site occupation number per site on a lattice), there 
appears an extinction threshold for the predator population that separates the 
absorbing state of a system filled with prey from the active coexistence regime 
through a continuous phase transition |47j . 

These examples all call for a systematic approach to include stochastic 
fluctuations in the mathematical description of interacting reaction-diffusion 
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systems that would be conducive to the apphcation of field-theoretic tools, and 
thus allow us to bring the powerful machinery of the dynamic renormalisation 
group to bear on these problems. In the following, we shall describe such 
a general method 021 149L 1301 which allows a representation of the classical 
master equation in terms of a coherent-state path integral and its subsequent 
analysis by means of the RG (for overviews, see Refs. [101 1111 IT^ ). 



2.2 Field theory representation of master equations 



The above interacting particle systems, when defined on a c?-dimensional lat- 
tice with sites i, are fully characterised by the set of occupation integer num- 
bers Ui = 0, 1, 2, . . . for each particle species. The master equation then de- 
scribes the temporal evolution of the configurational probability distribution 
P{{ni}; t) through a balance of gain and loss terms. For example, for the binary 
annihilation and coagulation reactions A + A —^ with rate A and A + A A 
with rate A', the master equation on a specific site i reads 



dPin,;t) 
dt 



A(n, -t-2) (n, + l)P(...,n, + 2,...;t) 

+A' {ni + l)n,P{...,n, + l,...;t) 
-{\ + \')n^ {n.,-l)P{...,n„...-t) , 



(154) 



with initially P({ni},0) = Y\iP{ni)^ e.g., a Poisson distribution P{ni) ~ 
tIq' e~"°/ni!. Since the reactions all change the site occupation numbers by 
integer values, a second- quantised Fock space representation is particularly 
useful [481 1491 1311] . To this end, we introduce the bosonic operator algebra 



ai,aj 



= = 



al,a\ 



ai, a 



(155) 



From these commutation relations one establishes in the standard manner 
that Qi and a| constitute lowering and raising ladder operators, from which 
we may construct the particle number eigenstates \ni), 



Oj \ni) = Hi \ni - 1) , a] \ni) = {n^ + 1) , aj a,; \ni) = \ni) 



(156) 



(Notice that we have chosen a different normalisation than in ordinary quan- 
tum mechanics.) A state with Ui particles on sites i is then obtained from the 
empty vaccum state |0), defined through |0) = 0, as the product state 



(157) 



To make contact with the time-dependent configuration probability, we 
introduce the formal state vector 



(158) 
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whereupon the Hnear time evolution according to the master equation is trans- 
lated into an 'imaginary-time ' Schrddinger equation 



dt 



The stochastic quasi- Hamiltonian (rather, the time evolution or Liouville op- 
erator) for the on-site reaction processes is a sum of local terms, -ffrcac = 
iJi(a|, Oi); e.g., for the binary annihilation and coagulation reactions, 

H,{ala,)^-X(l-af)aj-X'(l-aDalal . (160) 



The two contributions for each process may be physically interpreted as fol- 
lows: The first term corresponds to the actual process under consideration, 
and describes how many particles are annihilated and (re-)created in each 
reaction. The second term gives the 'order' of each reaction, i.e., the num- 
ber operator a] ai appears to the fcth power, but in normal-ordered form as 

a] af, for a fcth-order process. Note that the reaction Hamiltonians such as 
P60|l are non-Hermitean, reflecting the particle creations and destructions. In 
a similar manner, hopping between neighbouring sites (ij) is represented in 
this formalism through 

Our goal is of course to compute averages with respect to the configu- 
rational probability distribution P({ni}]t); this is achieved by means of the 
projection state {r\ = (Oinje"', which satisfies (VlO) = 1 and (P|a| = ("Pi, 
e"' Sij . For the desired statistical averages of observables that 



smce 



must be expressible in terms of the occupation numbers {ri.;}, we then obtain 
im) = E Fi{n.})Pi{n.ht) = {V\Fi{ala,}) |<|.(i)) . (162) 

Let first us explore the consequences of probability conservation, i.e., 1 = 
{V\^{t)) = {V\e-"*\<P{0)). This requires {V\H = 0; upon commuting e^."' 
with H, effectively the creation operators become shifted a| 1 + al, whence 
this condition is fulfilled provided Hi{a\ 1, ai) = 0, which is indeed satisfied 
by our explicit expressions (|160() and H161|l . By this prescription, we may also 
in averages replace a\ ai a^, i.e., the particle density becomes a{t) — (a^), 
and the two-point operator a^Oi a'^aj Oi Sij -f aj. 

In the bosonic operator representation above, we have assumed that there 
exist no restrictions on the particle occupation numbers n,; on each site. 
If, however, there is a maximum n,; < 2s -I- 1, one may instead employ a 
representation in terms of spin s operators. For example, particle exclusion 
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systems with = or 1 can thus be mapped onto non-Hermitean spin 
1/2 'quantum' systems. Specifically in one dimension, such representations in 
terms of integrable spin chains have proved a fruitful tool; for overviews, see 
Refs. [511 152L IS^ I54j . An alternative approach uses the bosonic theory, but 
encodes the site occupation restrictions through appropriate exponentials in 
the number operators e~°i°' |55| . 

We may now follow an established route in quantum many-particle theory 
| 56| and proceed towards a field theory representation through construct- 
ing the path integral equivalent to the 'Schrodinger' dynamics (|159l) based 
on coherent states, which are right eigenstates of the annihilation operator, 
o-i \4'i) = 01 with complex eigenvalues 0^. Explicitly, one finds 



ki) = exp 



4:«n |0) 



(163) 



satisfying the overlap and [over-) completeness relations 

d?6 



exp 



n 



(164) 

Upon splitting the temporal evolution (|159|l into infinitesimal steps, and in- 
serting Eqs. (|164|l at each time step, standard procedures (elaborated in detail 
in Ref. ^]) yield eventually 



{F{t)) (X / \{V[cl,m4>:]F{{4>i\)e 



(165) 



with the action 



Y^(^-^^(tf)+JJdt 



dt 



no<^*(0) , (166) 



where the first term originates from the projection state, and the last one from 
the initial Poisson distribution. Notice that in the Hamiltonian, the creation 
and annihilation operators a\ and are simply replaced with the complex 
numbers 0* and respectively. 

Taking the continuum limit, 4ii{t) ip{x,t), (j)*{t) ip{x,t), the 'bulk' 
part of the action becomes 



A[^i>, V'] = / d'^x / dt 



(167) 



where the hopping term (|161|l has naturally turned into a diffusion propaga- 
tor. We have thus arrived at a microscopic stochastic field theory for reaction- 
diffusion processes, with no assumptions whatsoever on the form of the (in- 
ternal) noise. This is a crucial ingredient for nonequilibrium dynamics, and 
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we may now use Eq. (|lt)7|) as a basis for systematic coarse-graining and the 
renormalisation group analysis. Returning to our example of pair annihilation 
and coagulation, the reaction part of the action (|157|l reads 

KcaclV-, ^) = -A (l - ^j^) V' - A' (l - V-) V- V'' , (168) 

see Eq. ()160|l . Let us have a look at the classical field equations, namely 
6 A/ Sip = 0, which is always solved by = 1, reflecting probability conserva- 
tion, and SA/S'tJj = 0, which, upon inserting ip = 1 gives here 

^^^^ ^DV i^ix, t) - (2A + A') i^ix, tf , (169) 

i.e., essentially the mean-field rate equation for the local particle density 
tjj{x,t), see Eq. ()149|l . supplemented with diffusion. The field theory action 
p67|l . derived from the master equation H154(l . then provides a means of in- 
cluding fluctuations in our analysis. 

Before we proceed with this program, it is instructive to perform a shift 
in the field ip about the mean-field solution, tp{x, t) ^ I + ip{x, t), whereupon 
the reaction Hamiltonian density (|168(l becomes 

Wrcac(V^, V-) = (2A + A') + (A + A') ^2 ^2 (^^q) 

In addition to the diffusion propagator, the annihilation and coagulation 
processes thus give identical three- and four-point vertices; aside from non- 
universal amplitudes, one should therefore obtain identical scaling behaviour 
for both binary reactions in the asymptotic long-time limit |57j . Lastly, we 
remark that if we interpret the action ^^[■0,?/'] as a response functional (|34|l . 
despite the fields ■0 not being purely imaginary, our field theory becomes for- 
mally equivalent to a 'Langevin' equation, wherein additive noise is added to 
Eq. (|169|l . albeit with negative correlator L[tp] — — (A + X')ip'^, which repre- 
sents 'imaginary' multiplicative noise. This Langevin description is thus not 
well-defined; however, one may render the noise correlator positive through 
a nonlinear Cole-Hopf transformation -ip — e'' , tjj — e~^p such that ipip — p, 
with Jacobian 1, but at the expense of 'diffusion noise' (x D p (Vp)^ in the ac- 
tion |58| . In summary, binary (and higher-order) annihilation and coagulation 
processes cannot be cast into a Langevin framework in any simple manner. 



2.3 Diffusion-limited single-species annihilation processes 

We begin by analysing diffusion-limited single-species annihilation k A 
I^EH]. The corresponding field theory action H167|) reads 



Aitp, V] = / d'^a; / dt 







d_ 

dt 



DV^ MA- A 



(171) 
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which for fc > 3 aUows no (obvious) equivalent Langevin description. Straight- 
forward power counting gives the scahng dimension for the annihilation rate, 
[A] = which suggests the upper critical dimension c?c(fc) = 2/(fc— 1). 

Thus we expect mean-field behaviour ^ (A see Eq. (|151ll . in any 
physical dimension for fc > 3, logarithmic corrections at dc = 1 for fc = 3 and 
at dc = 2 for fc = 2, and nonclassical power laws for pair annihilation only in 
one dimension. The field theory defined by the action H171|l has two vertices, 
the 'annihilation' sink with fc incoming lines only, and the 'scattering' vertex 
with fc incoming and fc outgoing lines. Neither allows for propagator renor- 
malisation, hence the model remains massless with exact scaling exponents 
77 = and z = 2, i.e., diffusive dynamics. 




Fig. 6. Vertex renormalisation for diffusion-limited binary annihilation A + A ^ 



In addition, the entire perturbation expansion for the renormalisation for 
the annihilation vertices is merely a geometric series of the one-loop diagram, 
see Fig.Elfor the pair annihilation case (fc = 2). If we define the renormalised 
effective coupling according to 

= ^ (172) 

where Bkd = kl r{2~d/dc) dc/k'^^^ (47r)''/'^<= , we obtain for the single nontrivial 
renormalisation constant 



Z-^ ^ 1 + -^^^^ (173) 



D {dc ~ d) 

to all orders. Consequently, the associated RG beta function becomes 

d_ 



.9i^.=-^(d-4+<?i?) , (174) 



with the Gaussian fixed point = stable for d > dc{k) = 2/(fc — 1), leading 
to the mean-field power laws (I151|l . whereas for d < dc{k) the flow approaches 

g* = 4(fc) - d . (175) 

Since the particle density has scaling dimension [a] — /i'*, we may write 
aij(/x, Dfl, no, 5_r) = /x'^ afl(Dfl, tiq /i^'', g/j), where we have retained the de- 
pendence on the initial density uq. Since the fields and the diffusion constant 
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do not renormalise (7d = and 7„o — —d), the RG equation for the density 
takes the form 

aB{no^i-'',gR) ^0 , (176) 

see Eq. (I75II . With the characteristics set equal to iii — {D t)^^/^, the solution 
of the RG equation 1176(1 near the IR-stable RG fixed point g* becomes 

aRino,t)^{DfiH)-''^'aR{rioiDfiH)^/^,g*) . (177) 

Under the RG, the first argument in Eq. (|177|1 flows to infinity. One therefore 
needs to establish that the result for the scaling function a is finite to all 
orders in the initial density 'S^, TT . One then finds the following asymptotic 
long-time behaviour for pair annihilation below the critical dimension 57, 55|. 

fc = 2 , d<2 : a{t) - (Dt)-'^^^ . (178) 

At the critical dimension, g{£) — > logarithmically slowly, and the process is 
still diffusion-limited; this gives 

fc = 2,d = 2: a{t) ^ (Dty^lniDt) , (179) 
fc = 3,d=l: a{t) ^ [(DtyHniDt)]^^'^ . (180) 

2.4 Segregation for multi-species pair annihilation 

In pair annihilation reactions of two distinct species A + B 0, where no 
reactions between the same species are allowed, a novel phenomenon emerges 
in sufficiently low dimensions d < d^, namely particle segregation in separate 
spatial domains, with the decay processes restricted to sharp reaction fronts 
on their boundaries [501 ■ Note that the reaction A + B ^ (d preserves the 
difference of particle numbers (even locally), i.e., there is a conservation law 
for c{t) = a{t) — b{t) = c(0) [SI]. The rate equations for the concentrations 

d{t) ^ -Xa{t)b{t) =b(t) (181) 

are for equal initial densities a(0) = b{0) solved by the single-species pair 
annihilation mean-field power law 

a(t) = 6(i) - (At)-\ (182) 

whereas for unequal initial densities c(0) = a(0) — 6(0) > 0, say, the majority 
species a{t) Coo = c(0) > as t — > oo, and the minority population 
disappears, b{t) 0. From Eqs. (|181|l we obtain for d > dc = 2 the exponential 
approach 

a{t)~a^^b{t)^e-'''-°'^^* . (183) 



""0 ^ — 



d 
dgn 
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Mapping the associated master equation onto a continuum field theory 
Hlt)7|) . the reaction term now reads (with the fields ■0 and Lp representing the 
A and B particles, respectively) [5^ 

7^,cac(?/',V','^,¥') . (184) 

As in the single-species case, there is no propagator renormalisation, and 
moreover the Feynman diagrams for the renormalised reaction vertex are of 
precisely the same form as for A + ^ ^ 0, see Fig. El Thus, for unequal 
initial densities, c(0) > 0, the mean-field power law ~ A i in the exponent of 
Eq. H183|l becomes again replaced with {Dt)'^/'^ in dimensions d < ~ 2, 
leading to stretched exponential time dependence, 

d<2: \nb{t) 1'^/^ d=2: \nb{t) t/ln{Dt) . (185) 

However, species segregation for equal initial densities, a(0) — 6(0), even 
supersedes the slowing down due to the reaction rate renormalisation. As 
confirmed by a thorough RG analysis, this effect can be captured within the 
classical field equations (SH- To this end, we add diffusion terms (with equal 
diffusivities) to the rate equations (|181l) for the now local particle densities, 

^ - DV^^ a{x,t) = -Xa{x,t)b{x,t) = -DV^^ b{x,t) . (186) 

The local concentration difference c{x, t) thus becomes a purely diffusive 
mode, dt c(a;, t) = D V'^c(a;, t), and we employ the diffusion Green function 

Goiq,t) = e(t)e"^'^^* , Go{x,t) = ^_M_e-V4Z.* , (137) 

compare Eq. (|23(l . to solve the initial value problem, 

c{x,t)= jd'^x Ga{x-x\t)c{x',Q) . (188) 

Let us furthermore assume a Poisson distribution for the initial density cor- 
relations (indicated by an overbar), a{x, 0) a{x', 0) = a(0)^ + a(0) 6{x — x') = 
b{x,Q)b{x',Q) and a{x,Q)b{x' {0) = a{Of , which implies c(a;, 0) c(cc', 0) = 
2 a(0) 6{x—x'). Averaging over the initial conditions then yields with Eq. 1] 188(1 



c{x,tY = 2a{Q) J d'^x'Goix-x^tf = 2a{0){8TTDt)-'^'^ ; (189) 

since the distribution for c will be a Gaussian, we thus obtain for the local 
density excess originating in a random initial fluctuation. 



a(0) 



\c{x,t)\ = \l-c{x,t)^ = 2 \j ^ {Sir Dty^^ . (190) 
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In dimensions d < — A these density fluctuations decay slower than the 
overall particle number ~ for d > 2 and ^ i"''/^ for d < 2 in a homo- 
geneous system. Species segregation into A- and _B-rich domains renders the 
particle distribution nonuniform, and the density decay is governed by the 
slow power law (PU^ . a{t) - b{t) - {Dty^/^. 

For very special initial states, however, the situation can be different. For 
example, consider hard-core particles (or A — > oo) regularly arranged in an 
alternating manner . . . ABABABABAB ... on a one-dimensional chain. The 
reactions A + B ^ % preserve this arrangement, whence the distinction be- 
tween A and B particles becomes meaningless, and one indeed recovers the 
power law from the single-species pair annihilation reaction. 

Let us at last generalise to q-species annihilation Ai + Aj — > 0, with 
1 5: * < J !i 9; with equal initial densities ai{0) as well as uniform diffu- 
sion and reaction rates. For g > 2, there exists no conservation law in the 
stochastic system, and one may argue, based on the study of fluctuations 
in the associated Fokker-Planck equation, that segregation happens only for 
d < ds{q) = 4:/ {q— 1) 63 . In any physical dimension d>2, one should there- 
fore see the same behaviour as for the single-species reaction A + A ^ 0; this 
is actually obvious for q = oo, since in this case the probability for particles 
of the same species to ever meet is zero, whence the species labeling becomes 
irrelevant. In one dimension, with its special topology, segregation does occur, 
and for generic initial conditions one flnds the the decay law 

a,(0~i-"(''^+Ct-i/2 , = (191) 

2q 

which recovers a{2) = 1/4 and a{oo) — 1/2. Once again, in special situations, 
e.g., the alignment . . . ABC D ABC D ABC D ... for g = 4, the single-species 
scaling ensues. There are also curious cyclic variants, for example if for four 
species we only allow the reactions A + B 0, B + C 0, C + D 
0, and D + A 0. We may then obviously identify A — C and B — D, 
which leads back to the case of two-species pair annihilation. Generally, within 
essentially mean-field theory one finds for cyclic multi-species annihilation 
processes ai{t) ^ t~"'-'^''^\ where for 

Remarkably, for five species this yields the borderline dimension 4(5) = 

for segregation to occur, hence nontrivial decay exponents a(5, 2) = i(-\/5— 1) 

in d = 2 and a(5, 3) = j{V5 — 1) in d = 3 that involve the golden ratio 

2.5 Active to absorbing state transitions and directed percolation 

Let us now return to the competing single-species reactions A — )■ (rate n), 
A ^ A + A (rate a), and, in order to limit the particle density in the active 
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phase, A + A ^ A (rate A). Adding diffusion to the rate equation (|152|l . we 
arrive at the Fisher-Kolmogorov equation of biology and ecology |46j , 

^f^^^^-D{r-V^)a{x,t)~Xa{x,t)\ (193) 

where r = {k — a)/ D. As discussed in Sec. 12.11 it predicts a continuous tran- 
sition from an active to an inactive, absorbing state to occur at r = 0. If we 
define the associated critical exponents in close analogy to equilibrium critical 
phenomena, see Sec. (|l.l|l . the partial differential equation H193|) yields the 
Gaussian exponent values 770 = 0, vq = 1/2, zq = 2, and = 1 = /3o- 

By the methods outlined in Sec. 12.21 we may construct the coherent-state 
path integral H167|l for the associated master equation, 



j A'^x Jdt ^j(^.^^-D V^^ V - (1 - V-) 

+a (l - ^ - \ (l - ifj iP ^l;^ 



(194) 



Upon shifting the field ip about its stationary value 1 and rescaling according 
to ip{x, t) = 1 + ^/(tJXS{x, t) and ip{x^ t) = ^/Xja S{x, t), the action becomes 



A\S, S] = A^x &t 



S [ ^ + D ir - V'') \ S - u [ S - S] S S + 



dt 

(195) 

Thus, the three-point vertices have been scaled to identical coupling strengths 
u — A, which represents the effective coupling of the perturbation ex- 
pansion, see Fig. |S1 below. Its scaling dimension is [u] = , whence we 
infer the upper critical dimension = 4. The four-point vertex oc A, with 
[A] = /i^"'', is thus irrelevant in the RG sense, and can be dropped for the 
computation of universal, asymptotic scaling properties. 

The action H195(l with A = is known as Reggeon field theory |f)5j . and 
its basic characteristic is its invariance under rapidity inversion S{x,t) <-!■ 
—S{x, —t). If we interpret Eq. (|195|l as a response functional, we see that it 
becomes formally equivalent to a stochastic process with multiplicative noise 
{{({x,t)) — 0) captured by the Langevin equation |66l I67| 

= -D(r- V^) Six, t)-u S{x, tf + C(a;, t) , (196) 
at ^ ' 

(C(a;, t) C(a;', t')) = 2u S(x, t) 8{x - x') 5{t - t') (197) 

(for a more accurate mapping procedure, see Ref. ) . which is essentially a 
noisy Fisher-Kolmogorov equation H193|l , with the noise correlator H197|l en- 
suring that the fluctuations indeed cease in the absorbing state where (5*) = 0. 
It has moreover been established ,69;, .7171 ITT] that the action H195|l describes 
the scaling properties of critical directed percolation (DP) clusters [TJj, illus- 
trated in Fig. Indeed, if the DP growth direction is labeled as 'time' t, 
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we see that the structure of the DP clusters emerges from the basic decay, 
branching, and coagulation reactions encoded in Eq. H194|) . 




Fig. 7. Directed percolation process (left) and critical DP cluster (right). 



In fact, the field theory action should govern the scaling properties of 
generic continuous nonequilibrium phase transitions from active to inactive, 
absorbing states, namely for an order parameter with Markovian stochastic 
dynamics that is decoupled from any other slow variable, and in the absence 
of quenched randomness [7S] . This DP conjecture follows from the follow- 
ing phenomenological approach to simple epidemic processes (SEP), or 
epidemics with recovery j46) : 

1. A 'susceptible' medium becomes locally 'infected', depending on the den- 
sity n of neighboring 'sick' individuals. The infected regions recover after 
a brief time interval. 

2. The state n = is absorbing. It describes the extinction of the 'disease'. 

3. The disease spreads out diffusively via the short-range infection 1. of neigh- 
boring susceptible regions. 

4. Microscopic fast degrees of freedom are incorporated as local noise or 
stochastic forces that respect statement 2., i.e., the noise alone cannot 
regenerate the disease. 

These ingredients are captured by the coarse-grained mesoscopic Langevin 
equation dtn — D (V^ — n + C, with a reaction functional R[n\, and 

s stochastic noise correlator of the form L[n] — nN[n]. Near the extinction 
threshold, we may expand R[n] = r + un -|- . . ., N[n] = v -\- . . and higher- 
order terms turn out to be irrelevant in the RG sense. Upon rescaling, we 
recover the Reggeon field theory action (|195|l for DP as the corresponding 
response functional 

We now proceed to an explicit evaluation of the DP critical exponents to 
one- loop order, closely following the recipes given in Sees. II. Hill. 51 The lowest- 
order fluctuation contribution to the two-point vertex function r^^'^\q,Lu) 
(propagator self-energy) is depicted in Fig.|H{a). The Feynman rules of Sec. II. 31 
yield the corresponding analytic expression 

ra^^)(g,.)^i. + i.(. + ,-) + gy^^ i./2i. + .l,V4 + fe- • ^'''^ 
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(a) 



(b) 



Fig. 8. DP renormalisation: one-loop diagrams for the vertex functions 
(propagator self-energy), and (b) _r(^'^) = — /^(^.i) (nonlinear vertices). 



a) r 



(1,1) 



The criticality condition /^'^'^^(0, 0) = at r = Tc provides us with the 
fluctuation-induced shift of the percolation threshold 



^ 1 



Oiu") . (199) 



D^Jk Tc + k 

Inserting t = r — into Eq. (|198|l . we then find to this order 

(g, ^) = i. + + g^) - !g / /-f^ + -- + f/4 , (200) 

and the diagram in Fig. ISl^b) for the three-point vertex functions, evaluated 
at zero external wavevectors and frequencies, gives 

r(-)({0}) = ^n^^m}) = -2. (l - ^ ^ ^-^) . (201) 



For the renormalisation factors, we use again the conventions (|56|l and 
but with 

UB^ = Z,uA'J%'-''-*y' . (202) 

Because of rapidity inversion invariance, Z~ — Zs- With Eq. (|5I-{|I the deriva- 
tives of _r(i^i) with respect to w, q^, and t, as well as the one-loop result for 
7^(1.2) jj^ Eq ||201() . all evaluated at the normalisation point tr = 1, provide 
us with the Z factors 



= 1 - 77^ ' ^" = 1 - 77^ • 203 

From these we infer the RG flow functions 

7s = vr/2 , 7c = -vr/A , 7r = -2 + ivn/A , (204) 
with the renormalised effective coupling 
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whose RG beta function is to this order 

Pv^VR[-e + SvR + Oivl)] . (206) 

For d > dc — 4:, the Gaussian fixed point = is stable, and we recover 
the mean-field critical exponents. For e = 4 — d > 0, we find the nontrivial 
IR-stable RG fixed point 

V* = e/S + Oie"^) . (207) 

Setting up and solving the RG equation Ht)5|) for the vertex function pro- 
ceeds just as in Sec. 11.51 With the identifications H83|) (with a = 0) we thus 
obtain the DP critical exponents to first order in e, 

r7 = -^ + 0(e2), l = 2-^ + 0(e2), ^ = 2-±+Oie'). (208) 

In the vicinity of v* , the solution of the RG equation for the order parameter 
reads, recalling that [S] = n'^^'^, 

{Sb,{tr, t)) « A*'/' ^^'^-^-^^/^ S (tr C'^ , v*R, Dr ^2 f+i'o t) , (209) 

which leads to the following scaling relations and explicit exponent values, 

/3 = ^ = l4 + 0(.^), a^^^l-Uoie^). (210) 
2 6 zv A 



Table 2. Comparison of the DP critical exponent values from Monte Carlo simula- 
tions with the results from the e expansion. 



Scaling exponent 


d = 


= 1 


d = 


-- 2 


d = 


-- 4 - 


e 






i 1.100 


V f= 


s 0.735 


= 


= 1/2 


+ e/16 + 0(e2) 


tc ^ t, ~ \T\ 


Z R 


i 1.576 


Z R 


:i 1.73 


2 = 


-- 2 - 


e/12 + 0(e2) 






s 0.2765 




s 0.584 


fi-- 


= 1 - 


e/6 + 0(e2) 




a ? 


s 0.160 


a f 


s 0.46 


a - 


= 1 - 


e/4 + 0(e2) 



The scaling exponents for critical directed percolation are known analyti- 
cally for a plethora of physical quantities (but the reader should beware that 
various different conventions are used in the literature); for the two- loop re- 
sults to order in the perturbative dimensional expansion, see Ref. [SHI- In 
Table 13 we compare the 0(e) values with the results from Monte Carlo com- 
puter simulations, which allow the DP critical exponents to be measured to 
high precision (for recent overviews on simulation results for DP and other 
absorbing state phase transitions, see Refs. |74ll75j V Yet unfortunately, there 
are to date hardly any real experiments that would confirm the DP conjec- 
ture |71[ I73| and actually measure the scaling exponents for this prominent 
nonequihbrium universality class. 
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2.6 Dynamic isotropic percolation and multi-species variants 

An interesting variant of active to absorbing state phase transitions emerges 
wiien we modify tiie SEP rules (1) and (2) in Sec l2.5l to 

1'. The susceptible medium becomes infected, depending on the densities n 
and m of sick individuals and the 'debris ', respectively. After a brief time 
interval, the sick individuals decay into immune debris, which ultimately 
stops the disease locally by exhausting the supply of susceptible regions. 

2'. The states with n ^ and any spatial distribution of m are absorbing, 
and describe the extinction of the disease. 

Here, the debris is given by the accumulated decay products, 

m{x,t) ^ K [ n{x,t')dt' . (211) 



After rescaling, this general epidemic process ( GEP) or epidemic with removal 
PS] is described in terms of the mesoscopic Langevin equation 



dS{x,t) 



-D (r -V^) S{x,t) ~ DuS{x,t) f S{x,t') Dt' + C(x,t) 

J —OQ 



dt 

(212) 

with noise correlator H197|) . The associated response functional reads [771 175] 

A[S,S]^ y'd'^xy'dt s(^-^^+D{r-V^)^S~uS^S + DuS J S{t') . 

(213) 

For the field theory thus defined, one may take the quasistatic limit by 
introducing the fields 

/oo 
S{x,t')dt' . (214) 
-oo 

For t — > oo, the action 12131) thus becomes 

^qst ['P,^]= / d'^x (filr - V'^ - u {>f - ip)] if , (215) 



which is known to describe the critical exponents of isotropic percolation [7^1 ■ 
An isotropic percolation cluster is shown in Fig. |^ to be contrasted with 
the anisotropic scaling evident in Fig. [^b). The upper critical dimension of 
isotropic percolation is dc = 6, and an explicit calculation, with the diagrams 
of Fig. [SJ but involving the static propagators Go{q) = l/{r + q^), yields the 
following critical exponents for isotropic percolation, to first order in e — 6 — d, 

ri = -^ + 0{e'), 1 = 2-^ + 0(6^), p ^ l ^ 1 + 0{e') . (216) 
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In order to calculate the dynamic critical exponent for this dynamic isotropic 
percolation (dIP) universality class, we must return to the full action (|213|l . 
Once again, with the diagrams of Fig. [SJ but now involving a temporally 
nonlocal three-point vertex, one then arrives at 

z^2-'-+0{^). (217) 
6 

For a variety of two-loop results, the reader is referred to Ref. [SHI- It is also 
possible to describe the crossover horn, isotropic to directed percolation within 
this field-theoretic framework |80[ 18 1| . 

Let us next consider multi-species variants of directed percolation pro- 
cesses, which can be obtained in the particle language by coupling the DP 
reactions Ai ^ (d, Ai ^ Ai + Ai via processes of the form Ai ^ Aj + Aj (with 
j i); or directly by the corresponding generalisation within the Langevin 
representation with {Q{x,t)) = 0, 

-g^^D, - RdS^]) S, + Q, R^[S^] =r,+J2 9v Sj + -- - , (218) 

j 

{Q{x, t)Q{x\ t')) = 2S,N,[S,] S{x - x') S{t - t') <5„- , N,[S.,] ^u, + ... (219) 

The ensuing renormalisation factors turn out to be precisely as for single- 
species DP, and consequently the generical critical behaviour even in such 
multi-species systems is governed by the DP universality class For ex- 
ample, the predator extinction threshold for the stochastic Lotka-Volterra 
system mentioned in Sec. 12.11 is characterised by the DP exponents as well 
|47|. But these reactions also generate Ai Aj, causing additional terms 
Y^j=^i9j^3 i'^ Eqs. H218|) . Asymptotically, the inter-species couphngs become 
unidirectional, which allows for the appearance of special multicritical points 
when several = simultaneously |82| . This leads to a hierarchy of order 
parameter exponents fik on the fcth level of a unidirectional cascade, with 

= /32 = i-^ + 0(e^),..., /3.. = ^-0(e); (220) 



52 Uwe Claus Tauber 



for the associated crossover exponent, one can show <l> — 1 to all orders |58 | . 
Quite analogous features emerge for multi-species dIP processes j58l 168) . 

2.7 Concluding remarks 

In these lecture notes, I have described how stochastic processes can be 
mapped onto field theory representations, starting either from a mesoscopic 
Langevin equation for the coarse-grained densities of the relevant order pa- 
rameter fields and conserved quantities, or from a more microsopic master 
equation for interacting particle systems. The dynamic renormalisation group 
method can then be employed to study and characterise the universal scaling 
behaviour near continuous phase transitions both in and far from thermal 
equilibrium, and for systems that generically display scale-invariant features. 
While the critical dynamics near equilibrium phase transitions has been thor- 
oughly investigated experimentally in the past three decades, regrettably such 
direct experimental verification of the by now considerable amount of theo- 
retical work on nonequilibrium systems is largely amiss. In this respect, ap- 
plications of the expertise gained in the nonequilibrium statistical mechanics 
of complex cooperative behaviour to biological systems might prove fruitful 
and constitutes a promising venture. One must bear in mind, however, that 
nonuniversal features are often crucial for the relevant questions in biology. 

In part, the lack of clearcut experimental evidence may be due to the 
fact that asymptotic universal properties are perhaps less prominent in ac- 
cessible nonequilibrium systems, owing to long crossover times. Yet fluctua- 
tions do tend to play a more important role in systems that are driven away 
from thermal equilibrium, and the concept of universality classes, despite the 
undoubtedly much increased richness in dynamical systems, should still be 
useful. For example, we have seen that the directed percolation universality 
class quite generically describes the critical properties of phase transitions 
from active to inactive, absorbing states, which abound in nature. The few 
exceptions to this rule either require the coupling to another conserved mode 
IHSIIMI; the presence, on a mesoscopic level, of additional symmetries that pre- 
clude the spontaneous decay A —> as in the so-called parity- conserving (PC) 
universality class, represented by branching and annihilating random walks 
A {n 1) A with n even, and A + A ^ [HHl (for recent developments 
based on nonperturbative RG approaches, see Ref. 'EE'); or the absence of any 
first-order reactions, as in the (by now rather notorious) pair contact process 
with diffusion (PCPD) 87 , which has so far eluded a successful field-theoretic 
treatment ISHi • A possible explanation for the fact that DP exponents have not 
been measured ubiquitously (yet) could be the instability towards quenched 
disorder in the reaction rates |89| . 

In reaction-diffusion systems, a complete classification of the scaling prop- 
erties in multi-species systems remains incomplete, aside from pair annihila- 
tion and DP-like processes, and still constitutes a quite formidable program 
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(for a recent overview over the present situation from a field-theoretic view- 
point, see Ref. ^2]). This is even more evident for nonequilibrium systems in 
general, even when maintained in driven steady states. Field-theoretic meth- 
ods and the dynamic renormalisation group represent powerful tools that I 
believe will continue to crucially complement exact solutions (usually of one- 
dimensional models), other approximative approaches, and computer simula- 
tions, in our quest to further elucidate the intriguing cooperative behaviour 
of strongly interacting and fluctuating many-particle systems. 
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